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ABSTRACT 

Emission intensity measurements of Hel and Hell lines have 
been made from a helium plasma in a longitudinal magnetic field 
(0-700 gauss) at neutral particle density 'v.lO 14 cm -3 . The lines 
include Hel transitions n 3 D+2 3 P (n=3, 4;>5 2, n 3 S-*2 3 P (n=3,4, 5), 
n 3 P-v2 3 S (n=3 / 4) / n 1 D-»-2 1 P (n=3, 4, 5,6), n 1 P^2 1 S (n=3,4), n 1 S-+-2 1 P 
(n=4,5) and Hell transition n=4 ->-n=3. Diagnostic measurements 
are made using Langmuir probes, triplet-singlet emission intensity 
ratios involving the transitions 4 3 S-»-2 3 P and 4 1 S-»-2 1 P / and applica- 
tion of population coefficient r ratios. As the magnetic field 
is varied in the 0-700 gauss range, the axial electron density 
changes from -ulO 10 cm" 3 to-vLO 12 cm" 3 . The electron temperature 
falls from vLl eV to ^6 eV. 

A collisional-radiative model is used for explaining the 
observed emission intensity enhancements in the magnetic field. 

The model consists of completely separated singlet and triplet 
sublevels of Hel up to n=5 and combined levels for n=6-12. Helium 
ionic levels up to n=13 are also incorporated. In addition to 
usual electron- imp act excitation-deexcitation processes, atom-atom 
collisions for Hel states and direct ionization to excited Helirni 
ionic states from all the Hel states are included. He X. and Hell, 
equations are solved simultaneously to yield populations of various 
states. The results show that the model explains the observed 
emission enhancements satisfactorily, and in this regard it is 
considerably better than some other models reported in the literal 
ture. Further, the model used in the present work also explains 
hel and Hell excited state populations in some recent helium 
plasma experiments reported by other workers. 



Chapter I 


INTRODUCTION 


In astrophysical plasmas, densities of particles in 
various excited states is determined by both physical and 

chemical processes. Since the excited ■.'States control-. - '-: ~ v' : 

line 

th e A emissivity of a plasma, measurement of line or continuum 
radiation often is a good probe for information about the 
plasma processes. A thorough quantitative understanding of 
these processes however is many times more difficult than a 
qualitative impression of the undergoing phenomena. Plasmas 
simulated under laboratory conditions often is the best 
recourse one finds to investigate plasmas of astrophysical 
interest. There is now a great deal of interest not only on 
the natural plasmas but also on plasmas to '.^realize fusion 
conditions and laser plasmas. While processes that are of 
major importance in fusion plasmas are not the plasma emissi-- 
vities, though the latter contribute to plasma cooling, but in 
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laser plasmas, afterglow plasmas, in all types of discharge 
plasmas, in beam-plasma and laser-plasma interactions detailed 
studies on population of excited states are of major interest. 

The present work is involved with emission of line radiation 
from a helium plasma and correlation of plasma conditions with 
excited state populations. 

Although it is true that emission of radiation from 
plasmas is perhaps the most easily noted plasma phenomenon, the 
emission processes do not necessarily constitute the most 
important processes among the plasma rate processes. What 
however is important, that for plasmas which are remote, as 
is true for astrophysical plasmas, emissivities are means by 
which not only the plasma constituents can be identified but 
are often the sole means by which the plasma conditions can 
be investigated. Had it been true that only a handful of pro- 
cesses are involved in excitation phenomena, it would have 
been indeed rather simple to understand plasma conditions and 
processes from the emissivities. But the processes which are 
responsible for the atomic excited state population at a given 
time are numerous a,nd their rates are dependent in a complex 
way on plasma conditions. The usual particles being the neutrals* 
positive ions, negative ions, and electrons, the excited state 
populations of heavy particles are determined by atom-atom, 
atom- ion, atom-electron collisions, ion-electron recombination 
processes, spontaneous emission processes, and, if the photon 
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flux in the plasma is significant, then by photo absorption 
processes as well. The population in an atomic bound state is 
thus controlled by a multitude of rate processes and though 
under particular conditions the relative importance of various 
processes change, even in the simplest of cases it is a formi- 
dable task to find an unambiguous answer about the sources of 
excitation and causes of population decay. In laboratory 
plasmas, the density and energy distribution of electrons in 
the plasma almost always play a dominant role in the excitation 
processes, particularly so when the photon flux is negligible. 
While theoretical emissivities can be calculated regardless of 
whether the actual energy distribution of plasma particles 
follows a simple Maxwellian pattern, categorization of plasmas 
on the basis of the degree and nature of equilibration attained 
is quite often used and this makes understanding of plasma 
excitation-deexcitation processes considerably simpler. 

If a plasma is contained within isothermal walls at 
temperature T and all its properties can be uniquely described 
by this temperature, then the plasma is said to be in complete 
thermodynamic equilibrium 1 (GTE). In GTE the distribution of 
energies of all the plasma constituents are uniquely derivable 
from the information on temperature and detailed information on 
the plagna system is not essential. From information on the 
temperature, determination of both the state of matter and 
radiation is uniquely possible. Under CTE conditions the laws 
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of Maxwell, Boltzmann and Saha are operative, from which energy- 
distribution of particles and population of particles in 
various states can be derived. The spectral distribution of 
plasma radiation under CTE condition is controlled by Planck's 
law. 

A more limited case of equilibration takes place when 
the criteria of local thermodynamic equilibrium (LIE) is 
satisfied. Under these conditions the electron velocity 
distribution is Maxwellian and this can be designated by an 
electron temperature T . In LTE, the plasma equilibrium is 
entirely determined by electron collisions and hence population 
of heavy particle energy states are compatible with T g satis- 
fying Boltzmann and Saha distributions, and the radiated 

energy is given by Kirchoff's law. For LTE to be satisfied, 

2 

a criterion in terns of electron density and temperature is 

n > 1.2 x 10 12 T ^ x(p,q) 3 cm -3 (l.l) 

e 0 

where x(p,q) is the excitation potential from level p to level 

q in eV. For example, when applied to the helium system at 

T =5 eV, at n > 10 ^^ LTE condition holds, 
e e 

When however the population of the energy states of heavy 
particles cannot be described by the T q values alone, one ha s 
non-LTE plasmas. For non- LTE plasmas, the population of 
excited states under a set of plasma conditions cannot be 
easily determined and one needs to consider all the excitation 
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and deexcitation processes in detail, and, depending on the 
physical conditions, transport processes like diffusion also 
need to be given due consideration. For non-LTE theoretical 
calculations, one therefore uses models in which collisional 

and radiative processes are taken into account in various 
degrees. 

In a typical such non-LTE model the processes that need 
to be taken into account, for calculations on excited state 
populations, are (i) electron Impact excitation, deexcitation 
and ionization, (ii) ion-electron recombination processes, i.e., 
continuum to bound state transitions either through three-body 
recombination or by radiative recombination, and (iii) photo- 
excitation, deexcitation, and ionization. Each of these 
processes has its characteristic energy dependent cross-section, 
which can be used to obtain a characteristic rate coefficient if 
the energy distribution of the colliding particles is known. 

As the Importance of the electron impact excitation-deexcitation 
processes depend on the electron density n and electron 
temperature values, the relative importance of photon Induced 
transitions in a plasma depends on the ambient photon flux and 
energy.. In this regard, one can thus further classify the 
plasmas as optically thick or optically thin. In optically thin 
plasmas, radiation does freely escape from the plasma and no 

reabsorption of the emitted radiation occurs. Optical escape 
3 , 

factors of all radiation can then be taken as unity. In 
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optically thick plasma, reabsorption of the emitted radiation 
takes place., and hence for such plasmas photo excitation must 
be an integral part of the non-LTE model. Here, the escape 
factors vary depending on the transition, the lowest value 
arises for the resonance lines of the system, and is numeri- 
cally taken as zero in the case of complete absorption. 

-I o 

At lower electron densities (n g ^10 ) a 'steady state 

2 

corona model 1 has been sometimes used to explain plasma line 
radiation from solar corona and also from low density labora- 
tory plasmas. In this model, which assumes optically thin 
condition, the balancing processes are collisional excitation 
ionization from the ground state on one side, and, the radia- 
tive recombination and spontaneous decay on the other. However, 
this simplified model is not very satisfactory because, even 
at electron densities 10 l0 -iO 13 , collisional excitation- 
deexcitations of excited states and three-body recombinations 
contribute significantly to the population densities. So a 
theoretical model which accommodates all the collisional and 
radiative processes taking place between various quantum 
states becomes desirable. Such a 'collisional-radiative model', 
which is of much wider applicability than the corona model, 
can be expected to give mere realistic population densities 
of excited states. 

The first detailed collisional-radiative model on both 
optically thin and optically thick hydrogenic plasmas was 
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^ 10VCX7J e.tb- 

proposed by^ Bates, Kingston and McWhirter. 4 The theory 
involves the following processes: 


N 

N 

N 

N 


( Z-l ) + + 

T <(p.^c), 
e 'K(c,pj 

H Z+ + e 

(z " 1)+ (p)+ 

K(p,q) ? 

e "k( q ,pf 


Z Hr ^ 

3(p) 

e ^ 

N (Z - l)+ (p) 

(z - i)+ ( P ) 

MP/qW 

N (Z - l)+ (q ) 


+ e 


+ e 


+ hv 


+ hv 


( 1 . 2 ) 

( 1 . 3 ) 

( 1 . 4 ) 

( 1 . 5 ) 


where K(p,c) is the ionization rate coefficient from the pth 
quantum state, K(c,p) is the rate coefficient for three-body 
recombination, K(p,q) is the rate coefficient of excitation 
(deexcitation) from the pth to the qth quantum state, B(p) 
is the rate coefficient for radiative recombination resulting 

in populating the pth state, A(p,q) is the spontaneous transi- 
tion probability from state p to state q. Neglecting transi- 
tions due to atom-atom, atom-ion or ion-ion collisions, 
assuming uniform distribution among degenerate levels, ignoring 
diffusion effects, for optically thin plasmas one can write 


h (p) 

n E TpT 


+ 


+ 


P(p) {n e [K(p,c) + l K(p,q)] + £ A(p,q>| 

W .q<P 


I p(q) n K(p,q) + T p(q) 
q^P q>P 


n E (q) 


A(q,p) + n e K(p,c) 


-2 - 

n^ 

_e 6(p) 

Xn E (p ) 


( 1 . 6 ) 
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where P(p) = n(p)/n.g(p), X = n e /n(N^ + ). rtg(p) is the number 
density of the species ^ + j_ n level p in Saha equilibrium 

with free electrons at density n g such that 

X n E (p)/n 2 = p 2 ( h 2 /2 w mlcT) 3 ^ exp (I At) (1..7) 

Ip being the ionization potential of the pth state. Eq. (l,6) 
should also contain a diffusion term if such effects are 
important. By substituting numerical values in eq. (1.7) 

Bates and others^ have shown that for a wide range of plasmas 

z+ 

« n g/ n(N ) for p^l, violated only in dense plasmas. 
However, if this condition is valid and the mean thermal energy 
is significantly less than the first excitation energy then 
f£ ( p) « n(l) v/hen the steady state condition is reached. In 
transient plasmas n(p) (p^l) values reach quasi-equilibrium 
values very rapidly after which the change in n(p) is insigni- 
ficant as the plasma decays. Under such condition fi(p)s e(p^l) 
can be equated to zero and a set of infinite linear equations 
can be obtained. For homogeneous steady state plasmas n(p)=0 
is directly obtained. The diffusion term goes to zero as there 

is no density gradient. For inhomogeneous stationary state 
5 

plasmas the diffusion term has to be incorporated if diffusion 
loss of n(p) is significant. For solution of eqij..{li,^5 , an 
infinite matrix is avoided by grouping together all the levels 
with high quantum numbers above p and designating them 3 d y a 
such that p(o) is close to unity. The latter is true, for 
higher states where the collisional processes are relatively 
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more Important than the radiative ones, and n(p) almost satis- 
fies Saha equilibrium densities. If p(l) is considered as 
known, a solution of the set of equations (1.6) can be 
expressed as P(p) = r Q (p) + r^p) p(l), r Q (p) and r^p) being 
the population coefficients of the pth state and are functions 
of only n^ and T q . This means that from a knowledge of r Q (p) 
and r^(p) the population of the pth quantum state can be 

determined if n and T are known. 

0 0 

Determination of r Q (p) and r^(p) requires information on 
atomic parameters like rate coefficients, transition probabi- 
lities etc. Therefore, the accuracy to which the population 
coefficients can be obtained depends on the accuracy to which 
the atomic parameters are known or can be calculated. For 
hydrogenic plasmas, McWhirter and Hearn^ worked out the 
formalism for calculation of r^Cp), r^(p) coefficients using 
reduced variables 6= T/Z 2 , n(c) = n e /Z 7 and n(p) = Xn(p) /Z 11 
On hydrogen plasmas, Johnson and Hinnov have recently made 
extensive calculations based on a new set of electron colli- 
sional cross-sections and computed the values of r Q (p) and 
r 1 (p) in the electron density range 10 5 -10 18 and electron 
temperature range 250°-8.2 x 10 6 °K for the states p=2, 3, 4, 5, 6. 

Whereas the collisional -radiative model is simplest to 
apply in the case of atomic hydrogen plasmas, the need for 
quantitative understanding is far more general, which means 
application to non-hydro genic plasmas. Clearly the complexities 



10 


are enormous when the energy levels multiply many times and 
various sublevels are to be considered, which implies that 
drastic simplifications in the model are essential in order 
to make the problem tractable. Before discussing the work 
hitherto done on the problem of helium plasmas, a brief 
mention about some of the recent applications of the colli- 
sional radiative model to non-hydrogenic plasmas will be 

informative. 

3 

Park has calculated spectral line intensities in a non- 
equilibrium nitrogen plasma using the method of Bates, Kingston 
and McWhirter.^ The objective of his work has been to derive 
populations of atomic nitrogen states as a function of 
electron density, electron temperature, and then use intensi- 
ties of any two atomic nitrogen lines to determine the electron 
temperature of any nitrogen plasma. In his model, he grouped 
some of the bound states of nitrogen atom, and a total of 41 
groups were involved. He showed that if the equilibrium 
relation is used for temperature determination, non-equilibrium 
effects cause the apparent to be higher than the true in 
an expanding plasma. Above 8000° K, in an expanding non-equili- 
brium plasma, the departure from the equilibrium line intensity 
ratio becomes negligibly small. 

9 

Recently Giannaris and Incropera have studied a cylin- 
drically confined argon arc plasma and applied the collisional- 
radiative model to obtain the excited state populations. Here 
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also, the precise energy level structure of Argon was not 

retained and levels were grouped, according, to common electron 

of prin cipal quan tuETniiaib e r ® 

configuration, and level s^ higher than 15” were assumed to be 

in equilibrium with free electrons in the plasma. Assuming 

that the net absorption and emission rates equal for resonant 

transitions, and using known n , n and T from thermochemical 

e o e 

properties of Argon arc, local values of the rate coefficients 
and radial distribution of the level populations were determined. 
They found that for the most severe condition of the study, at 
the wall of 35 A arc, excitation non- equilibrium effects extend 
to the 6f level. For all arc regions other than the near-wall 
region, non-equilibrium effects extend only to the 4p level. 
Further, the non- equilibrium effect is always an overpopulation 
of the levels relative to the Boltzmann-Saha equilibrium value. 

Helium plasmas are interesting for several reasons. 

Relative to hydrogen, helium is the first atom where the atomic 
sublevels play an important part in plasma excited state popula- 
tion distribution, yet helium being, the simplest of the multi- 
electronic atoms, the system is relatively most convenient to 
handle. At the same time, -the upper' quantum states in helium 
resemble the atomic hydrogen states more closely than any other 
atomic system. It is also the simplest atomic system with 
metastable states (2^3 and 2 s) which do not have optical 
transitions to the ground state and hence requires a collisional 
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mechanism for depopulation to the lowest state. Helium is 
also known to show the widest deviation from LTE under a set 
of electron density and electron temperature conditions. 
Experimentally, helium being mono-atomic, it is a relatively 
simple matter to ascertain its neutral concentration in plasmas, 
whereas hydrogen being a diatomic gas, it is often difficult 
to determine the percent atomic concentration in hydrogen 
plasmas where the gas is only partially dissociated. Further, 
collisional-radiative model already worked out for hydrogenic 
ions can be easily applied to He + ions in a helium plasma. 

A mention of some of the helium plasma experiments will be 

appropriate here though all of than are not directly related to 

the present work. Most of the experiments in which helium 

excited state population have been the subject, of study used 

afterglows. In Princeton plasma physics laboratories Motley 
10 

and Kuckes studied He afterglow discharges in the B-l stel- 
larator and showed, using microwave phase shift measurements 
for determination of n^ and conductivity measurements for T , 
that the principal mechanism of charge removal in low tempera- 
ture highly ionized helium plasma is the three-body recombina- 
tion reaction. Hinnov and Hirechberg ' measured absolute 
intensities from the spectrum of afterglow obtained from B-l 
stellarator and their conclusions are compatible with those of 
Motley and Kuckes. Further, Hinnov and Hirschberg derived the 
recombination coefficients as a function of n Q in the range 
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1115 

10 -10 and in the T g range 0.1-0. 4 eV. They also studied 

the population densities, both theoretically and experimentally, 
from the afterglow emission intensities, and found that the 
experimental population densities compare well with the theory, 
taking only n=3,4,5 levels. Their calculations, however, did 
not take into account splitting of sublevels. 

•I O 

Johnson has done experiments on afterglow and also on 

12 13 

ohmic heating of discharges in the n g range of 10 -10 and 

T of 0.25 and 4-14 eV. He used a collisional-radiative model 
e 

with levels up to n=25 with separated sublevels up to n=8. 

He assumed that above n=25 Saha equilibrium is valid. He 
also calculated collisional-radiative recombination coeffi- 
cient which he compared with experimental values. Johnson and 
Hinnov 1 ^ have done experiments on afterglow with ohmic heating 
in the density range 10 12 - 10 12 and T g range 0.04-1 eV. Comparison 
with experimental results shows that the rate coefficient 

K(2 5 S *► l^s) cannot adequately account for the observed popula- 

3 1 

tions, a more likely path is through the process 2 S -*-2 S 
followed by deexcitation to the ground state. It may be 
mentioned here that in an earlier work Myers 15 had suggested 
that the mechanism of metastable cooling is through atom-atom 
collisions, ionization, followed by ion-molecule reactions. 

In Collins and Robertson's 16 flowing afterglow experiments, 

at n of 5 x 10 12 and maximum temperature of 0.15 eV, it was 
e 

found that above n=6 Saha equilibrium exists and below this 
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level the population densities can be explained by the colli- 

sional-radiative model. They also concluded that the primary 

P r " oCess *-)f populating these levels is the collisional— radiative 

recombination of He ions. Gusinow, Gerardo and Verdeyen 1 " 1 

studied afterglow plasma at 2-8 Torr in the n e range lO 13 — 10^~* 

and T g range 0.09 to 0.6 eV and used He— Ne laser interferro— 

metric technique to measure n and T . They found that the 

experimental value of a is comparable to that calculated from 

the collisional— radiative model of Bates, Kingston and McWhirter. 

Born and Buser studied afterglow phenomena from a condenser 

discharge in the pressure range 0.3-0. 7 Torr, n g = 10 12 -10 14 , 

T e = eV. According to them the decay of the metastable 

3 

state 2 S was not by usual deexcitations; diffusion to the 
wall is a more important cause for the observed temperature 
decay. 

Of other kinds of experimental systems, Robben, Kunkel 

and Talbot 19 have done experiments in helium plasma jet at 

T 0 range 0.15-0.3 eV and n £ range of 10 12 -10 14 and has shown 

that K(p,q), (p=q+l ) , values remain constant in the range of 

experimental conditions used and the measurements are consis- 

20 

tent with Gryzinski's values. Va Fugol, Pakhomer and Reznikov, 
using high frequency discharge experiments at 1 Torr, T g 12-4 eV, 
have shown that the metastable 2 3 S lifetime is 2x1 0“ 4 sec and 
the metastables are responsible for He* formation as well as 
for the intense afterglow. Newton and saston 22 conducted 


21 
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experiments on pulse discharges in the n g range 8x1 O^-l. 8x10^ 
and pressure range Oi82-l*2 Torr. Their major conclusion is 
that the metastable atoms deactivate by volume deexcitation 
in the n e range 10 14 -l 0 15 , but at lower densities wall deactiva- 
tion is likely. 

Ikee and Takeyama 23 have recently studied brush cathode 
plasmas in a magnetic field. They determined n g and T g from 
measurements of spectral intensities of the continuum spectrum 
which follows the series 2 3 S+-n 3 P. From a plot of log (e(X) /° c ) 
vs 1/ X(e is the emission coefficient and a c is the photo- 
ionization cross-section) they determined T q/ and from the 
absolute value of e(X) the value of n . At a discharge current 
of 500 mA and a pressure of 0.9 Torr, T g is 0.17 eV and n g 
1 . 8x1 0 1 3 . The recombination rate calculated from the experi- 
mental results, which is equal to 2.2x10 10 cm 3 sec 1 at 

T = 0.17 eV and n * 1.8 x 10 13 , is in good agreement with 

e e 

that calculated from the collisional-radiative recombination 
theory of Bates, Kingston and McWhirter. The n Q was found to 
’ increase by one order of magnitude compared to n Q without the 
magnetic field. 

I^ee and Takeyama 24 have also studied population densities 

of Hel excited states in brush cathode plasmas. They measured 

11 11 

spectral line intensities for the transitions 2 Sen P, 2 Pen D, 
2 3 sen 3 P, 2 3 Pen 3 S and 2 3 P«-n 3 D up to n=10 and determined the 
population densities of the excited states. They found that 
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for n>,5 the states are in thermal equilibrium and from Boltzmann 
plot determined the temperature, which was found compatible with 
the free electron temperature measured from recombination conti- 
nuum. By interpolation of Saha increment values b as given by 
Drawin and Emard / theoretical population densities were 
calculated for n^, values of the experiment. A comparison 
of the theoretical population densities with those obtained from 
the experiments on absolute intensity measurements show that 
the results match within the errors of interpolation of b^ and 
the experimental measurements. 

Recently Otsuka, Ikee and Ishii 26 have made studies on 
helium plasma using their TPD (Test Plasma by Direct current 
discharge) machine. They made spectroscopic measurements of 
HeX and Hell line intensities at various points of the magnetized 
plasma column and also measured plasma parameters at the points 
of emission intensity measurements. The T g range of their 
experiments was 0.2-8 eV, n g ^10^^ and system pressure 'vlO - 2 
Torr. From emission intensities, absolute population densities 
of the He I states and Hell (n=4) states were determined. 

Estimates were made of n ++ densities and these were applied to 
predict, using a collisional radiative model, the expected 
population densities for Hell (n=4) state. They also proposed 
a mechanism, for the observed phenomena in TPD machine, accord- 

, -f- f- 

xng to which He recombination is the major process upstream 
in the plasma column causing Hell excited states to be populated 
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by cascading of higher states. At downstream conditions, with 
lower electron temperatures, recombination of He + ions consti- 
tutes the major process causing higher levels of neutral helium 
to increase in population density which results in enhanced 
emission from such states. 

Emission enhancement studies of Hell line 4686 S. have been 

27 

made by Takeyama and Takezaki in a longitudinal magnetic 

field in the pressure range 0.4-4 Torr and magnetic field 

intensity of 0-6 kilogauss. They find that the emission 

intensity has a maximum enhancement dependent on the pressure. 

The characteristic pressure for maximum intensity increases as 

the longitudinal field increases. They, however, do not offer 

any explanation for the emission enhancement. In a later work 

28 

Takeyama and Takezaki measured enhancement of helium neutral 

i 

lines and found that the relative intensities belonging to 

2 1 P«-n 1 D (6678, 4922, 4388 A) and 2 1 S«-n 1 P (3965 and 3614 A) 

series of Hel increase as the longitudinal field is increased. 

The enhancement becomes maximum at about 2 kilogauss. The 

11 3 3 

lines belonging to the other four series 2 P<-n S, 2 S^n P, 1 
3 3 3 3 

2 P«-n S and 2 P*-n D (except the line 5876 A) however do not ' 

show any enhancement, on the contrary, the relative intensities 
monotonically decrease as the field increases. 

While some of the experimental work reported above involved 

calculations, there has been considerable work solely from 

29 

theoretical side. D'Angelo calculated the value of recombination 
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coefficient a in heliun systems in the range 0.26-0.9 eV and 
11 13 

n e range 10 -10 using a model of three body recombination 

and found the calculated values matched well with the large 

values of a obtained in Motley and Kuckes^ experiments. Bates 

and Kingston used approximate expressions for He -e elastic 

scattering cross-section to calculate quasi-equilibrivm 

electron, ion temperatures and the corresponding collisional 

radiative recombination coefficient a and compared the 

calculated values, in the n range 10 9 - 10 15 , with the experi- 

12 

mental results of Hinnov and Hirschberg on afterglows. Bates, 

31 

Bell and Kingston have also calculated the properties of a 

decaying optically thick plasma in the T g range 0.09-1.4 eV 

8 16 

and n e range 10 -10 in which they took into consideration 
Penning ionization, and concluded that the primary process 
which determines the population of the first excited state is 
collisional-radiative recombination. 

Drawin and coworkers have made extensive calculations of 

32 

population densities in helium systems. Initially, Drawin 

made some calculations of population densities of He I states 

5 

and in a recent work Drawin, Emard and Katsonis have repeated 
the same calculations with refined atomic constants like term 
values, transition probabilities and cross-sections. The 
results are applicable to homogeneous stationary, transient 
and inhomogeneous diffusion dominated plasmas. The general 

■ 4 

model is similar to that of Bates, Kingston and McVJhirter but 
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for actual calculations they have taken a total of 51 levels/ 

1 11 
which include the l S ground state, and the states 2 S, 2 P, 

3 3 

2 S, 2 P. The remaining 46 levels are equally divided between 
singlet and triplet series without special distinction between 
the sublevels, S> P, D etc. of the same principal quantum 
number, and going upto n=25. In the earlier work atom-atom 
exchange collisions between singlet and triplet levels had been 
incorporated, but in the recent work they are not included; - 
which helps keep the rate equations linear. Instead, they have 
added electron-atom exchange collisions in the model and this 
partially makes up for the error involved in neglecting the 
atom-atom exchange interactions and, for degrees of ionization 

_3 

larger than 10 , the electron-atom collisions are much more 

efficient than atom-atom collisions. They made calculations 
for the limiting cases of optically thin and optically thick 
plasmas and also for intermediate cases where the degree of 
radiation trapping is intermediate between the extreme cases. 

As indicated earlier, in homogeneous plasmas the rate 
equations simplify easily. Treatment of inhomogeneous or 
transient plasmas can be simplified if (a.) for transient plasmas 
the relaxation times are short, and (b) for the inhomogeneous 
plasmas, the relaxation lengths are smaller than the smallest 
lateral dimension of the plasma in the direction of the gradients. 
In case of the latter, the system can be taken to be in quasi- 
homogeneous state and the divergence of the diffusion fluxes of 
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all the excited levels can be put equal to zero. Drawin and 

others'’ have given results of population density calculations 

in the electron density range lO^-lO 2 ^ and T^ range 0.09 to 

11 eV for optically thin, optically thick and partially 

optically thick helium plasmas using explicit values of 

escape factors for various transitions involved. They also 

provide results for Saha decrements of ground state population 

and electron densities, and they find that trapping of resonance 

photons make significant changes in the ground state Saha 

decrements for electron densities <10 . In all the cases 

of radiation trapping the Saha decrements are greater than 

17 

unity for n g <10 

Drawin and others have also applied the results of their 

calculations mentioned above to the experiment of Boersch, Geiger 

33 

and Topschowsky on He arc discharge. The latter experimenters 

did not provide a clear value for the electron temperature. 

Using the collisional-radiative model, Drawin and others^ made 

calculations, for a range of T g values and for a ground state 

population of 5.5x10 , for the population of excited states 

using both inhomogeneous and the homogeneous stationary state 

model. In the n range 10 12 -10 13 and T = 2.8 eV good agreement 

0 €2 

between the measured and calculated values were found. 

The above sample list of experiments and theoretical calcu- 
lations on the helium systsn shows that work in both the cate- 
gories has been numerous and widely varied, with overlaps in some 



21 


cases. Part of the difficulty in comparing one set of results 
with another however is that, in spite of special care, it is 
difficult to completely specify all experimental conditions in 
discharges, and sometimes unspecified parameters contribute 
significantly to the observed phenomena. But in principle, 
all experimental results should be explicable from one physical 
model, when appropriate allowances are made for different 
experimental conditions. No experiment or theoretical work 
referred above has explicitly undertaken such an extensive 
task. This is probably because a great deal more of experi- 
mental data is necessary, which implies that the understanding 
at this stage has necessarily got to be evolutionary in 
character, and gradually the gaps will be bridged. The present 
work was •undertaken as an effort in this direction. 

The purpose of the present study is to examine, experi- 
mentally, the applicability of the collisional-radiative model 
to low- temperature moderate density heliun plasmas. The way 
it differs from' some of the experiments reported earlier is 

that by varying some experimental parameters, a wide range of 
n , T conditions are realized in the plasma. The extended 
range of n , T provides a wider scope of testing this model. 
Specifically, the emission of He I and Hell line radiation is 
studied under varying n , T conditions, and attempts made to 
understand the observed line intensities from application of a 
collisional radiative model. 
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Experimentally, the system that is used in the present 
work is a stationary state plasma system which has similarities 
to the positive colunn. Whereas in general the plasma electron 
density can be increased by injecting more power and more 
electrons into the plasma column, one alternative and conve-r- 
nient way to effect this is to gradually inhibit diffusion of 
the electrons to the wall containing the plasma by using a 
confining magnetic field. If the plasma is ignited in a 
cylindrical column and a longitudinal magnetic field is 
applied, the effect on the plasma charged particles is to 
confine the latter preferentially near the axial region with 
a consequent increase in electron density. Thus a much larger 
and efficient increase of electron density can be achieved in 
this manner concurrent with the reduction of transverse 
diffusion. Though the quantitative aspects of such inhibition 
of transverse diffusion 34 3 ^ is a subject of considerable 
discussion and the exact dependence of the diffusion coeffi- 
cient of the magnetic field is not entirely clear, the 

monotonic increase of the density is an experimentally observed 

3 8 : 

fact. Also known is the fact that with increase in field, 

the electron temperature decreases due to more frequent 

collisions with electrons and consequent loss of electron 
38 

energy. Thus in a plasma column subject to a longitudinal 
magnetic field, with change in magnetic field intensity the 
plasma column behaves as a reaction volume which moves through 
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a coupled variation of electron temperature and electron 
density and thus controls the excited state populations through 
detailed excitation-deexcitation processes. While the diffusion 
processes are unavoidably present, in a stationary state plasma 
soon a radial concentration profile is established for all the 
plasma species. For electrons and other charged particles one 
expects a maximum concentration on the axis. From a theoretical 
analysis of diffusion processes in cylindrical plasma columns 
it follows that the electron density radial profile is like 

“D Q 

a zero order Bessel function R(r) = R(0)j (2.40& -) where 

o a 

r is the distance from the axis and a is the tube radius 
(2.2 cm in the present case). This means that for a consi- 
derably wide zone near the axis, the electron density can be 
taken as practically invariant, e.g., within 10% variation, 
if r = 0.35 cm (size of the baffle in the present ' experiments ) . In 
this region the plasma resembles a homogeneous stationary state 
plasma. That the excited state populations follow somewhat 

similar distributions has been experimentally shown by 

q 

Giannaris and Incropera. 

In the present work amission intensities of He I and Hell 
lines are measured end-on from the near-axial region of a 
plasma column subjected to a longitudinal magnetic field. 

Plasma diagnostic measurements are carried out also in the 
same region. With increase in the magnetic field, the line 
emission intensities increase, and these emission enhancement 
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factors are experimentally measured. These emission enhance- 
ment results are then compared with the predictions of a 
co 11 isional -radiative model/ for both optically thin and 
optically thick plasma conditions. For application of the 
collisional-radiative model to ' stationary state 
plasmas and particularly to the emission enhancement results 
of the present investigation/ independent calculations of 
population coefficients are also carried out. The results of 
these calculations are compared with those of other workers. 
Further/ the model used in the present work is applied to 
some helium plasma experiments reported in recent literature. 



Chapter II 
experimental 


All the experiments are carried out in a long cylindrical 
discharge' tube with conducting wall. The discharge is sustained 
by low energy electrons injected into the plasma column from 
an electron gun. The longitudinal magnetic field applied on 
the plasma is provided by a solenoid aligned coaxially with 
the discharge tube axis. The radiation from the plasma is 
taken out through a quartz window, collimated by a set of 
baffles external to the discharge tube and then analyzed using 
a grating monochromator. 

The details of the discharge tube are schematically shown 
in Fig. 2.1. The tube is about 1 m in length and 2.2 cm in 
internal diameter. At both ends of the tube, chambers of 
internal diameter 9 cm and length 10 cm are connected and are 
used for mounting the electron gun, optical windows, vacuum gauges. 
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pumping line connections, gas inlet valve and electrical feed- 
throughs. The central tube as. well as the side chambers, made 
out of 304 non-magnetic stainless steel, are provided with 
water cooling jackets. The system is pumped by a 10 cm oil 
diffusion pump backed by a mechanical pump. There is an 
isolation valve between the diffusion pump and the discharge 
tube which permits, during an experiment/ partial or total 
isolation of the system from the pumping line. 

The outer cover of the electron gun consists of a 6 cm 
long 2.1 cm diameter stainless steel tube. The electron gun 
assembly is mounted on a 5 cmi diameter stainless steel flange 
which itself is mounted on a 7.5 cm diameter backing plate 
using a set of 2.5 cm long ceramic insulators. A 0.5 mm 
diameter tungsten filament, shaped in the form of a 1 cm 
diameter circular ring is mounted on two 2 mm diameter tungsten 
rods, through which filament power is supplied. 

The filament is located 2 mm inside the top end of the ' '■ 
outer cover and inside a guard ring made of 1 mm tungsten wire 
wound in the form of a coil and mounted coaxially with the 
electron gun cover. The guard ring is electrically insulated 
from the outer wall of the gun and the filament. The top end 
of the gun has also a grid made of 0.25 mm tungsten wire placed 
2 mm away from the plane of the filament. The, grid which is 
in contact with the body of the gun, is also shorted to the 
discharge tube. The electron gun can be pushed inside the 
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discharge tube in such a manner that the backing plate rests 
on the electrical feedthroughs mounted on the collar of the 
central tube. There is a clearance of 1 cm diameter through 
the filament ring and the backing plate for end-on measurements 
of plasma radiation. 

The magnetic field is provided by a solenoid mounted 
coaxially with the central discharge tube. The solenoid is 
wound with 20 SWG super-enamelled copper wire on a 95 cm long 
4.2 cm diameter spool and is capable of generating about 
250 gauss per ampere. The maximum current capability of the 
wire, with ambient air cooling, is about 4 amp. 

The filament is biased 150 volts negative with respect 
to the discharge tube, the latter being kept at ground poten- 
tial. The guard ring is at the same potential as that of the 
filament to minimize electron emission towards the wall. At 
the other end of the discharge tube a 1.5 cm diameter electrode 
or a stainless steel wire gauge disc, kept at ground potential, 
was placed, which constrained the 90 cm long plasma column to 
be free of external field. No observable' difference in experi- 
mental results was found in these two modes of operation. » 

The electrical arrangement is shown in Fig. 2.2. The 
filament requires about 25 amp current (Ip) for 4 amp or 
electron emission current (i Q ) and this is provided by a 
variable dc power supply. The drop across the filament is 
negligible compared to 150 V bias voltage. The bias voltage 
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is provided by a 200 volt 5 amp floating power supply through 
a 0.5 ohm limiting resistor and an ammeter for measurement of 
the electron emission current injected into the plasma column. 
The filament is located entirely inside the magnetic field and 
away from the rest of the mounting assembly. This assists the 
emitted electrons to be drawn into the plasma column. A 300 
volts power supply is used for the solenoid and is connected 
through an ammeter and a 1 ohm (20 watt) resistor. The voltage 
sensed across this resistor provides a measure of the magnetic 
field applied. 

The discharge system is pumped by a 10 cm oil diffusion 
pump’ backedpby t mechhni-cdl -punp to a'ppesbure. of ; -10— mm of Hg-V- 
From a careful study of leakage rates it was concluded that 
lower pressures could not be reached mainly because of the 
inefficiencies of the pumps and not due to any significant 
leakage. The pressure is measured by thermocouple and Bayard- 
Alpert ionization gauges. Both the gauges and their controls 
were calibrated against a Stokes McLeod gauge with helium gas 
in the system. Helium gas (Matheson, specified purity 99.99%) 
is introduced into the system through a needle valve. The 

required gas pressure (all experiments reported in the present 

—3 

work were carried out at 5 x 10 - Torr) is maintained by 
controlling the flow rate through the valve. By concurrent 
control of the needle and isolation valves study could be done 
in a partially closed system to see the influence of flow rate 
on enhancements of radiation in magnetic field. 
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The radiation is observed along the axis of the discharge 
tube through the backing plate of the electron gun. To see 
that only parallel radiation from the plasma column along the 
axis is collected, the radiation is passed through two baffles 
of 0.7 cm diameter opening centred along the axis of the tube. 
The baffles are located such that they reject all radiation 
leaving the tube with an angle more than 0.5° with the axis. 
Measurements were also carried out with baffles .35 cm 
diameter without noticeable changes in experimental results/ 
except the absolute intensities. The collimated beam is then 
focussed on the slit of the monochromator by a concave mirror 
which is placed 1 focal length away from the slit. A high 
resolution spectrometer built in our labo rate ry 40 sis used as the 
monochromator (Appendix Fig. A. 1 ) . It uses a Bausch & Lomb 
123x154 mm replica plane grating with 1200 grooves/mm blazed 
at 7500°A and employs Ebert-Fastie mounting in Czemy-Tumer 
modification. The spectrometer can be operated either in 
2.1 m or 6. 4 m mode leading to a first order dispersion of . 

3.4 A/mm and 1.12 S/mm respectively. In the first order, the 
range of the spectrometer is 5040 to 12020 St. The spectrometer 
is equipped with a camera for use with photographic plates as 
well as a set of two photomultiplier detectors for simultaneous 
measurements of 2 line intensities, if the lines are closely 
located. The slit width used for most measurements was 
0.035 mm. The spectrometer was operated in 2.1 m mode for 
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the present experiments. The two photomultipliers (1P28, 
maximum response in 3400+500 A region) are mounted in two 
different housings, which have slit openings of .020-. 030 mm 
x 25 mm located on the focal plane of the monochromator 
focussing elements. The dynodes of the photomultipliers are 
biased from a highly stabilised Fluke power supply (0.25% 
stability) and normally operated at 1000 volts. The photo- 
multiplier anode currents are taken to an instrumentation 
amplifier operated in the difference mode as shown in Fig. 2.2. 
The amplified signal is taken to a digital voltmeter (sensi- 
tivity O.i mV) for direct readout or to an X-Y recorder. 

Enhancement Experiments 

The emission enhancement experiments were carried out as 

follows. The discharge system was first evacuated 'and 

helium admitted to the system through a needle valve to attain 

the required pressure, -150 volts bias was applied to the 

filament and the filament current was slowly increased till 

the desired emission current i was reached. The monochromator 

0 

was set for the required emission line such that the slit of 
one of the photomultipliers is on the centre of the line while 
position of the other is adjusted on the continuum about 50 A 
away from the emission line. The output of the difference 
amplifier is recorded (which subtracts for the continuum 
intensity). The magnet current is then sys ton atic ally increased 
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which permits measurement of line intensity at various magnetic 
field values. An intensity measurement in the absence of the 
field provides I q . Enhancement ratios I B /I Q can then easily 
be calculated. 

Emission Intensity Ratios at Two Emission Currents* (I /I ) 

i l i 2 B 

These measurements are similar to those for simple emission 
enhancements in magnetic field. Here the information desired is 
the ratio of line emission intensities at two emission currents, 
at various magnetic fields. The magnetic field is first set to 

a desired value; the emission intensity of a given line is 
measured first at one emission current and then at the other. 
Independent zero corrections for the intensities are made by 
moving the photomultiplier to continuum, because the continuum 
intensity changes significantly with change in the emission 
current. Due to limitations of the biasing power supply an 
increase in the emission current results in lowering of the 
bias voltage, so every such change required readjustment of 
the bias voltage to -150 volts. 

Langmuir Probe Measurements 

A 2 mm long single probe was used, made of 1 mm diameter 
tungsten wire. The probe was provided with a glass shield at 
the base of the wire to reduce effects due to sputtering 
deposits. The conditions (i) the dimension of the probe < mean 
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free path of the ions, (ii) probe current « discharge current, 
and (iii) negligible disturbance of the plasma field by the 
probe were satisfied. The probe was introduced into the column 
from the pump end and the probe tip was located at the centre 
of the column. The nature of variation of probe characteris- 
tics as the probe position was varied along the axis was noted. 

Measurement of He I and Hell Line Intensity Ratios 

In order to get an idea about the excited state popula- 
tions of Hell using the line intensity ratio method, HeT 4713 A 
and Hell 4635 A lines were used. For this pair of lines, with 
a fixed setting of the grating (2nd order, 2.1 m operation) 
the detectors were moved for measuring the intensities, one 
after another. Appropriate continuum corrections were also 
made. Similarly, for determination of electron temperature, 
intensity ratios of sing let- trip let lines 5047 A to 4713 A and 
4921 A to 4471 A were measured at various emission currents at 
different magnetic field value. 
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RESULTS 


Emission Enhancements in Magnetic Field (Ig/I^ 

n 3/ 1 L-* 2 3 ' 1 (L+l) series of He I lines were used for line 
intensity enhancement measurements in the presence of magnetic 
field. The details of the lines chosen are given in Fig. 3.1. 
In all, 16 Hel atomic lines were studied. They lie in tine 
range 3000-7000 A and are easily accessible in the monochro- 
mator used. For the series n 1 D-»-2 1 P (n = 3, 4, 5, 6), plots of 

I /I i.e., the ratio of the line intensity at a given 
B / o’ 

magnetic field B to that at zero magnetic field as a function 
of the field are shown in Figs. 3.2, 3.3, 3.4 and 3.5. For 
each line I B /I Q was measured at emission currents of 0.5 amp, 

1 amp, 2 amp and 3 amp. In Figs. 3.7, 3.8 and 3.9, variations 
of i /i with magnetic field are shown for the triplet series 
n 3 S+2 3 P- (n = 3,4,5). In Figs. 3.6 and 3.10, the (I B /I 0 -B> 




1 

Fsg.3.1 He lines chosen for measurement (energy 
levels not to scale , wavelengths are in X) 
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11^3 

plots for the series n D 2 P and n S **• 2 P are compared at a 
constant emission current of 1 amp. 

From the results shown in Figs. 3. 2-3. 5 and 3. 7 -3. 9, the 
following general observations can be made. 

(a) Lower the emission current, higher is the enhancement 
factor I B /l . 

(b) For the first line of any series, I B /I Q increases rather 
rapidly in the lower magnetic field range (0-350 gauss), 
and then the increase is relatively less steep. For the 
other lines of the series, the enhancement factor increases, 
passes through a rather broad maximum, and slowly decreases 
at the highest magnetic fields. 

(c) For a given line, as the emission current increases, the 
magnetic field at which a maximum in the enhancement 
factor is observed shifts to decreasing field values. 

This trend is somewhat less noticeable for higher lines 
of any series. 

(d) From Figs. 3.6 and 3.10, it can be seen that at a given 
emission current, the magnetic field at which the maximum 
enhancement is observed decreases with increase in the 
energy of the emitting state. Sharp fall in the enhance- 
ment factor is observed for higher lines in the series. 

Similar enhancement studies were carried out for all the 

n 3 D +2 3 P (n = 3 , 4,5), n 3 P -»-2 3 S (n = 3,4), 


four other series, viz.. 
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n S -*-2 P (n = 4/5), n^P-*-2^S (n = 3,4). For brevity, values 
°f ^B^o'^max' i*s*/ the maximum value of Ig/I attained, for 
all the lines of the above series and the corresponding 
magnetic field, at the four emission currents, are presented 
in Table 3.1. The Table also gives (I B /I ) ?00 , i.e., I B A 0 
value at 700 gauss. The general observations made above can 
also be verified from the Table. 

The only measurable Hell line from the plasma, at 4686 % 
(n=4-*-n=3), was used for emission enhancanent measurements, and 
the results are shown in Fig. 3.11. 

'.A '■■Figures'- 3.2 to 3.11 represent data from a single 
run. The range of variation of experimental I b /I q for the Hel 
lines (each measured at least 4 times , at time-intervals some- 
times of months) is given in Table 3.1. It is estimated that 
Ig/I values are dependendable within ± 1 %. 

Error in l_/I measured for 0.5 amps emission current is 
higher (^15%) than those at higher amission currents because 
of extremely small signal observed at zero field, in general, i^ 
■values - are dependable within +3%. Estimated variations in 
pressure are in the range 4 — 6x1 0 Tor r , but even when the 
pressure was increased up to 10 microns, changes in the Ig/I^ 
factors were found to he within the experimental error reported. 

At lower pressures, experiments could not be conducted 
because of the limitations of the filament power supply and 
also due to inadequate cooling arrangements of the discharge tube 
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On crlo-sing the isolation valve partially for a short time, no 
noticeable difference was found in the enhanconent factors. 
This shows that enhancement factors are quite independent of 
the gas flow rate through the tube. Number of experiments at 
i @ =3 amp was limited, as at high currents heating of the 
plasma tube becomes significant. 

i, i, 

Results of Ratio Measurements I„ x /I 0 

A typical plot of i.e. , the ratio of the line 

intensity at 3 amp emission current to that at 1 amp emission 
current at a given magnetic field, versus the magnetic field 
is shown in Fig. 3.12 for the series n D 2 P (n = 3, 4, 5, 6). 

In Fig. 3.13, the ratio of the line intensity at 2 amps to 
0.5 amp emission current for the same series is plotted as a 
function of the magnetic field. 

Similar measurements were made for other series and a 
typical plot for a triplet series is shown in Fig. 3.14 for 
a ratio of 3 amp to 1 amp emission current. These results 
are to be used later for electron density determination. Here 
again,- a few generalisations can be made: 

(a) The ratio is lower for higher lines at a given magnetic 
field. 

(b) Ratio for the same line decreases with increase in the 
magnetic field. 








m 


3 amp ; 1 amp 


3 3 P-*2 3 $ 

O— «• V 


-*> 3889 


S' ""Too' 
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400 
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“>■ 
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Fig. 3. 14 

Line intensity 

ratios 

at 3 amp 

to 1 amp 


emission current for 

n 3 P~ 2 3 S 

series 




51 


In Table 3.2, the line intensity ratios at two different 
emission current ratios (4 and 3) at zero and 700 gauss are 
given for all the lines studied. The observations made above 
on the basis of n D -*-2 P series generally hold good, as can 
be seen from Table 3.2. Studies were made with different 
sets of emission currents keeping the i ratios seme. When 
2.25 amp to 0.75 emission current ratio was chosen, the 
observed line intensity ratios were higher when compared with 
the values corresponding to 3 amp to 1 amp, at all magnetic 
field values. 

The ratio measurements were made in a separate set of 
experiments and the data given in Table 3.2 are reproducible 
within 6%. 


Temperature Measurements from Probe: 


* 


In Fig. 3.15, typical I-V characteristics of the probe 

at 1 amp of emission current are given. In Fig. 3.16 the plots 

\ 

of logarithm of electron current (In i ) against the probe 
voltage at various magnetic field and at 1 amp amission current 
is shown. Electron temperature obtained from the slopes of 
such curves at various magnetic field values is plotted in 
Fig. 3. 19. 

At zero magnetic field the probe current was extremely 
low and as such no dependable measurements could be made. At 
50 gauss, the typical characteristics of a single probe were 



i.i. S'. hAfKFUe 

cewikal LIBRARY 

— 
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TABLE 3.2 



Emission 

intensity 

ratios at 

two Sfiissibn 

currents 





(A/i 12 ) 

v B 7 B 




No. 

Wavelength Transition 
A 

il7i 2 =2 amp/0. 5 amp=4 
B=0 gauss B=700 guass 

i 1 /i_=2.25 
1 , 1 0.75 
Bat)' gauss 

amp/. 

amp=3 

B=7 OOgauss 

1 

5875 

3 3 

4.0 

2.60 

2.9 

2.10 

2 

4471 

3 3 

4 D+2 P 

3.6 

2.23 

2.5 

2.10 

3 

4026 

5 J D+2 3 P 

3.1 

2.02 

2.25 

1.5 

4 

5015 

1 1 

3 P+2 S 

3.0 

2.2 

2.2 

1.95 

5 

3964 

1 1 
4 x P-*-2 x S 

3.1 

2.18 

2.33 

1.82 

6 

5047 

4 1 S+-2 1 P 

4.0 

2.5 



n 

/ 

4437' 

5 1 S-+2 1 P 

4.0 

2.0 








i./i 2 = 3 amp/1 amp 






= 3 


8 

6678 

1 1 

3 D+2 P 

4.0 

2.62 

3.0 

2.32 

9 

4922 

1 1 

4 D+2 P 

3.96 

2.4 

2.93 

2.18 

10 

43 88 

1 1 

5 1 D-*2" L P 

3.85 

2.0 

2.80 

1.80 

11 

4144 

1 1 
6 X D*2 P 

3.8 

1.9 

2.72 

1.60 

12 

3889 

3 3 

3 p-^2 P 

4.0 

2.9 

3.0 

2.10 

13 

3187 

3 3 

4 P+2 S 

4.0 

2.86 

2.87 

1.96 

14 

7065 

3 3 

3 S+2 P 

4.2 

2.36 

3.00 

1.76 

15 

4713 

3 3 

4 S*2 P 

4.0 

2.48 

3.00 

1.45 

16 

4121 

3 3 

5 J S+2 P 

3.0 

2.55 ' 

2.35 

1.48 
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observed, i.e., ion as well as electron currents saturated. 

At higher magnetic field, the electron current did not reach 

saturation though ion saturation current could be measured. 

These observations are similar to those of some earlier 
38 

workers. At higher emission currents, the electron collec- 
tion current was high and probe heating took place. In this 
region also measurements could not be made. 

With the observations that could be made, the varia- 
tions of electron temperature were from 10.0 eV to 5.6 eV in 
the range of 50-700 gauss. On extrapolation to zero field, 

T° obtained was 11.4 eV, 
e 

A set of measurements was made at 0.5 amp emission 
current and the variation in the electron temperature in the 
range of 50 to 700 gauss was just about the same as shown in 
the T versus B plot in Fig. 3.19. 

Density Measurements from Probe Data 

Electron density was calculated using the equation 

J + = 0.4 n + e a (- M - e ) 1/2 (3.1) 

37 

given by Bobm and others, where J + is the ion saturation 
current in amperes, e is the electron charge in coulombs, 

M mass of helium ion, a • area of the probe and 3c is the 

4 * 

Boltzmann constant. n + (=n^is calculated using obtained 
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as above, variation of n g was found to be from 6.8x10^ cm - ^ 
to 4.5 x 10^^ over the magnetic field range 50 to 7-00 gauss 
as shown in Pig. 3.19. 

It is to be mentioned here that axial variation of T 

e 

was found to be within 10% at all magnetic field values. 

Further, at high and low fields, the ion saturation current 
was found to be practically invariant of the axial position. 

At medium magnetic field strengths the variation was found to 
be ^15%. with respect to the value at the centre of the colunn 
- with high values at the electron gun end and low values at 
the pump end. 

Electron Temperature from Emission Intensity Ratio : 

In Figs. 3.17 and 3.18 the experimentally observed 
intensity ratios of 5047 8- to 4713 and 4921 to 4471 % lines 'are 
plotted against magnetic field at different emission currents. 
It shows that the intensity ratios and therefore the electron 
temperature, is independent of emission current at a given 
magnetic field. The plasma electron temperature from 
Fig. 3.17 is obtained by the following method. 

It has been noticed that the electron impact excitations 
from the ground state to 4^S and 4^S states are strongly depen- 
dent on electron temperature. Considering the effects due to 

1 3 

cascading from the higher to 4 S and 4 S states and excitation 
to these levels from the ground state, the intensity ratio of 
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5047 S and 4713 & lines given by Sbvie 41 as 

*5047 = K _(1 1 S^ t.4 1 S) A (_4*S + 2 1 T?)_ 4-7/3 ' 

4713 K (1 1 S ^4 3 S) A (4 3 S +2 3 P)*S 7047 (3 * 2 ^ 

where h s are collislonal excitation coefficients, and 
A(p-^q) = A(p*q)/ \ A(p+ q ) where A’s are spontaneous • •• ■ 
transition probabilities. Values of K's and A‘ s used here 
will be described in detail in Chapter 4. Emission intensity 
tios calculated from eq. 3.2 at various electron temperatures 
is then plotted against T Q . A comparison of this plot and, 
the experimentally observed line intensity ratios yield T 
values which are shown in Fig. 3.19. (4 1 D- > -2 1 p)/(4 3 d^2 3 p) 

intensity ratios of Fig. 3.18 were not used for T q measure- 
ments because these are known to be far less dependable than 
the (4 S-*-2 P)/(4 3 S-»-2 3 P) line ratios. 

Effect of variation of Aperture Size on I /I 

B / o ’ 

In Fig. 3.20, a plot of I B /I o vs B for 3 889 St line at 
1 amp emission current is shown with different aperture sizes. 
As mentioned earlier baffle sizes varying from 0.3 cm to 0.7- cm 
diameter were used. To ensure collection of radiation from 
the axial region, 0.3 cm baffles were used at both ends of the 
tube for cfine al.-jgnment^ - which give an angular divergence 
of«0.5°* But ever, then the same results were obtained. 






Chapter IV 


ANALYSIS AND DISCUSSION 


The primary purpose of the present work was to study 
emission enhancements in helium plasmas in presence of a 
magnetic field and correlate the experimental results with 
the fundamental population-depopulation processes occurring 
in the plasma. in the previous Chapter results of the emis- 
sion enhancement studies for a number of Hel and Hell lines 
have been presented. In terms of the plasma parameters 
obtained, conservative estimates using eq. (l.l), or eq. ( 4 ) 
of Ref. 42, show that the plasma is under non-LTE condition 
throughout the entire range of magnetic field and emission 
current conditions of the experiment. Since explanation of 
observed emission intensities requires knowledge about 
population densities of excited states, and for non-LTE 
conditions theoretical population densities cannot be easily 
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obtained, a detailed collisional-radiative model applicable 
to inhomogeneous stationary state plasmas is necessary. 

In this Chapter, (i) a collisional-radiative model is 
formulated which takes into account series of levels from 
both He I and Hell systems; (ii) using this model, Hel and 
Hell population coefficient calculations are carried out for 
a wide range of n g , T g conditions, and the results presented; 
(iii) these results are then applied to the heliim plasma 
investigated in the present work and comparison made between 
the experimental emission enhancement factors and those 
predicted theoretically. Later, (iv) the present experi- 
mental results are reviewed in the perspective of the popu- 
lation coefficients calculated by Drawin and the rate 
coefficients of Park, and (v) the collisional radiative model 
used in the present work applied to Stellarator C helium 
plasma experiments of Johnson and TPD machine experiments of 
Otsuka, Ikee and Ishii. 

The Collisional Radiative Model Used in the Present Work: 

The basic collisional radiative model has already been 
described in Chapter 1. The details of calculations differ 
on. whether the system is optically thick or optically thin, 
the number of quantum levels talcen into account for detailed 
balancing, and the choice of values of the rate coefficients 
A(p,q), K(p,q), K(p,c), 8(p). For the Hel system, detailed 
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calculations for the excited state populations have been 
done by Johnson'*' 0 and . also by Drawin and Board. ^ Calcula- 
tions for the Hell system have been done by McWhirter and 
Hearn^ and also by Drawin and Emard.^ The Hel calculations 
of Johnson involve separated sublevels up to n=8 and upper 

limit of n=25 but their population densities are available 

43 

in limited T and n combinations. Drawin and Board's 
e e 

calculations are for a wide range of temperature and density 
conditions and also for different values of optical thickness. 
But they have used grouped singlet and triplet levels above 
n=2 . Experimental part of the present work uses Hel lines 
which involve individual sublevels of singlet and triplet 
states of n = 3, 4, 5, 6 and therefore none of the above 
calculations are directly usable. As regards Hell excited 
states, simultaneous ionization-excitation processes needed 
to be included in the model to explain the observed popula- 
tion- densities at the T , n range of the present experi— 
ments. Resort was therefore taken to an independent 
calculation of the non-LTE population densities of Hel and 
Hell excited states taking separated sublevels of Hel 
as -are' necessary for analyzing the emission enhancement data. 

The levels used for detailed balancing are shown in 
Table 4.1. The salient features of the scheme are that it 
uses separated sublevels up to n=5, and all the sublevels 
above n=5 are combined together. The highest level taken 



62 


TABLE 4.1 

Energy levels of He I used in the collisional-radiative model 


^ — w ‘ 

* - - v - 

- - - 

- - — - - - 


' - - - k r - - 


mvi * / 

Level No . 
(p) 

State 

Energy 
(cm~i ) 

... E P.„„ 

g 

P 

Level 

<P> 

No . State 

Energy 

; e p... . 


1 

l^S 

0 

1 

17 

3 1 

4 J F, 4 F 

191452 

• 28 

2 

2 3 S 

159856 

3 

18 

4 1 P 

191493 

3 

3 

2 1 S 

166278 

1 

19 

5 3 S 

193347 

3 

4 

2 3 P 

169078 

9 

20 

5 1 S 

193663 

1 

5 

2 1 P 

171135 

3 

21 

5 3 P 

193801 

9 

6 

3 3 S 

183237 

3 

22 

5 3 D 

193917 

15 

7 

3 1 S 

184865 

1 

23 

5 2 D 

193919 

5 

8 

3 3 P 

185656 

9 

24 

3 1 

5 J F / 5 1 F, 

5 3 G,5 1 G 

193921 

64 

9 

3 3 D 

186102 

15 

25 

S 1 ? 

193943 

3 

10 

1 

3D 

186105 

5 

26 

6 (n) 

195251 

144 

11 

3 1 P 

186210 

3 

27 

7 

196070 

392 

12 

4 3 S 

190298 

3 

28 

8 

196595 

512 

13 

4 1 S 

1909 40 

1 

29 

9 

196954 

648 

14 

4 3 P 

191217 

9 

30 

10 

197213 

800 

15 

4 3 D 

191445 

15 

31 

11 

197398 

968 

16 

4 1 D 

191447 

5 

32 

■ 12 

197543 

1152 


Energy 

levels 

of Hell 

used’ in 

the model 




i = 1 to 13 

E, = 109677 x 4 (l - ij) + 198311 cm' 



and 
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into consideration for He! is n=12. it will be shown ^ 

that with further incorporation of higher levels the calcula- 
ted population densities do not vary significantly. i„ 


addition. Hell levels 
total of 45 levels, 
thermally distributed, 


up to n=13 are included to give a 
Sublevels of Hell are considered to be 


The following processes are included in the model. 

The Hel states are denoted by p ' s and q ' s and Hen states by 
i ' s and j * s 


He (p ) + e <: K iP^g^> 

K(q,p) 

He (q) 

+ e 

(4.1) 

He (p) + e KCp^i) ^ 

He + (i) 

+ e + e 

( 4 . 2 ) 

K(i,p) . 

He (p ) + e < 

He* 1 * (l ) 

+ e + e 

(4.3) 

He(p) . A (p,q) 

P > q 

He (q) 

+ hv 

(4.4) 

He + (l)+ e 

He (p) 

+ hv 

(4.5) 

He(p) + HeClX^h 

. He (q) 

+ He (l ) 

(4.6) 

He (p) + He(l) < KN / (p / i)i=1 > 
KN (i, p ) i=l 

He + (1) 

+ He ( 1 ) + e 

(4.7) 

He + (i)+ e ) .> 

j # i T 

He + (j) 

+ e 

(4.8) 

He + (i)+ e 

I<(c, x) 

He ++ 

+ e + e 

(4.9) 
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He + (i) He + (j) + hv (4.10) 

He++ + e '^‘V> He + (i) + hv (4.1l) 

Most of the above processes are commonly used in colli 
sional-radiative models, Eq. (4.2) represents simultaneous 
ionization-excitation to Hell states and eqs. (4.6) and (4.7) 
describe neutral-neutral collisions for excitation and ioniza 
tion respectively. 

Since, in the T g range of the present study, contribu- 
tion even from the state HeII(l) towards the population 
density of Hel excited states is insignificant, contribution 
from Hell excited states to Hel levels is neglected. 

In the above scheme only spontaneous decay processes 
are included. For optically thick conditions, reduced 
transition probabilities A(p,q), A(i, j) are used to take 
into account reabsorption of radiation from plasma. This is 
equivalent to inclusion of reverse processes of eqs. (4.4) 
and (4.10). 

* 

Putting n/n + (l) = X and n /n ++ = Y, one has the 
0 0 

following equilibrium relations: 


K(p, q) .n E (p) = K(q,p) 

.n E (q) 

(4.12) 

(n /X).K(i,p) = K(p,i) 

n E (p) 

(4.13) 

i=l 


K(i,j) nj(i) = K(j,i) 

4 ( 5 ) 

(4.14) 

(n e /Y)K(c,i) = K(i,c) 

n E ( i ) . 

(4.15) 
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n(p) and n^(i) denote nunber densities of Hel and Hell states 

and n e , n ++ denote mnber densities of electron and doubly 

rged ion s . i^Cp) and n£(i) describe Saha equilibrium 

population densities respectively of Hel and Hell states 

given by eg. (1.7). Here c denotes the second continuum. 

For the above processes, time derivatives of excited state 

populations after incorporating equilibrium densities can be 
written as 


fi(p) 

n E (p) 


and 

ft + (i) 

n E (i) 


= - »<P> (n e [*f K(p,i) + 3 f K (p, q )] 


32 


q^p 


r 

+ n(l) I KN(p,q) + J A 

q^P p>q J 


32 


+ plq P(q) KN(p ' c J ) * n(l) ' n ^(p) + p(q)K( P/ q).n e 

+ I P(q).A(q,p) + n ,K(p, i) ,, , , , 

q>P n g lp; e . (i=i) + xr^Cp) ^(p/ 


+ no) K n C l, P) 

x nJLw 


(4.16) 


!+ <i) {n e [K(i,o) + 'f K(i,j)J + l A(1#j) , 


i^j 


13 


+ J P + (j> K(i,j).„ e + J + (J) “E 


nitjJ 


&3 

32 

I 

p=l 


J >1 


n+(i) 


* A Cx / x ) 


n F (p) 


J p(p ) . n . ' Kfn.i) 4- yi i r(A x-. N , 0 


n 


fid) 


s ‘ n+<l) K<P ' 1) + n - K(i ' C) + Yn+U) ■ 


[p(p) = n(p)/n E ip), » + (i)=n + (i)/n+(i) ] 


(4.17) 
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It is assumed that loss due to diffusion to the walls 
for all the excited states of Kel and Hell is negligible. In 
that case §qs. (4.16) and (4.17) represent the controlling 
equations for all Hel and Hell excited states. The eqs. (4.16) 
and (4.17), after equating the right hand side to zero 
(for p^l, i/4l), can be combined to give a set of 43 linear 
simultaneous equations. The inhomogeneous stationary state 
solution, in terms of the local ground state population 
densities of Hel and Hell (which have contributions from 
diffusion losses), can be written as 

p(p) = r Q (p) + r^(p) . p(l) (4.18) 

p + (i) = r*(i) + r*(i) p(l) + r 2 (i).P + (l) (4.19) 

, 9 n + / f ) 

It is to be noted that n E (p)/n^(i) = . exp ( (E -E .. )/kT Q ) , 

J i n P J - 

and in order to solve the set of linear equations, an explicit 

value of tSafes ratio n'(l)/n +l is required. The solution is 

therefore obtained in terms of n , T , T and n + (l)/n ++ , and 

S O Cf 

given these parameters, population densities of all the 
excited states can be determined in terms of n(l) and n + (l). 

In eq. (4.18), r (p) gives the contribution towards p(p) from 
the first continuum and r^(p). p(l) is the contribution from 
the ground state of Hel. In (4.19), r*(i) is the 

contribution from the second continuum towards p + (i), 
r^(i)p(l) gives the contribution from the ground state of 
Hel while r^ (i) < 3+ (l) gives the contribution from the ground 
state of Hell. 
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Effect of auto ionization processes is assumed to be 
relatively less significant. Population of He + ground state 
due to autoionization will not change the results of the 
present calculations since Hell (i ) and Hel(p) population 
densities ' are calculated on the basis of assumed values of 
He + ground state population (i.e., n + (l) = n e >. Further, 

radiative decay from the autoionizing states populate the 
very high excited states of He I preferentially. Atom-atom 
collisions are included in the general- programme although 
at neutral densities and electron temperatures of the present 
experiments it makes little contribution. 


Diffusive losses of the metastable states are neglected 

for the following reasons. 45 The total relaxation frequency 

for the metastable states can be written as v =1/T=n(l ) o^<v> 

+ K n(l)+n lK(m,q) where t is the relaxation time, n(l) is 
m eq ^ 

the ground state population density, a ^ is the momentum trans- 
fer cross-section, <v> the average velocity of the ground 


state particles, the rate coefficient for neutral -metastable 
collisions and K(m, q) the rate coefficient for excitation- 
deexcitation by electronic collisions. For worst case calcu- 


lations, at neutral densities of the present experiments, 

i.e., 2 x 10 14 , T 'ulOOO K, o, = 2.4 x 10“ 15 cm 2 , n o = 10 10 

g u e 

T e <vll eV, the second term of the above equation is negligible 

4 

and the first and third terms are of the same order, *v5xlO 
sec -4 . At all higher n g values, the third, term becomes so 
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dominant that diffusive losses of the metastable population 
density may be neglected in comparison with volvme deactiva- 
tion. It may be added here that, in some recent calculations 4 ^" 
also, metastable diffusive losses have been neglected. 


Choice of Atomic Parameters : 

The various collisional and radiative rate coefficients 
used in the calculations are obtained as follows. 


Spontaneous transition probability A(p,q) values for 
Hel are taken from the tables of Wiese, Smith and Glennon. 47 
A(i,j) values for Hell are obtained from the tables of Green, 
Rush and Chandler. 4 ^ 


For calculation of collisional ionization rate coeffi- 
cxents K(p,i; , Gryzinski's formulation is used. Ioniza- 
tion cross-sections for Hel states is then given by 

He + •>, 

He e Q <vy = E -2 “ «p<VV' (4 - 20) 

P 


where E^ is the ionization potential of the pth level and E g 
is the energy of the colliding electron. Putting E^/E = x, 
g^(x) is given by 


g n (x} = x ( x+l~ )3/2 l 1 + I (1 " ln (2 .?+(x-1)%(4.21) 


2x 


Ionization rate constant K(p, i) fair* the pth level is obtained 
by averaging the corresponding cross-sections over Maxwellian 
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distribution of velocity of the colliding electron, i.e.. 


K(p, i) 


r" He 

/ He e^ A A /E ) f(v).v.dv, 

v* e P 


(4.22) 


where 

f (v)dv = (2/ir)^ (m/kT e ) 3/2 exp(-% mv 2 /kT e ) . v 2 .dv 

and v* = (2E /m)A With the simplification h mv 2 = E , and 
P e 

substitution of the values of the atomic constants. 


K(p, i) 


2.918 
T 3/2 


/ g (x) e ^ x dx, 
1 p 


where a = E p AT g . This can also be written substituting for c^(x), 

K (p, i ) = (2.918/T e 3/2 ) exp (-a) f“f (z) e" Z dz (4.23) 

o 

after putting (x-l) = y and ay = z. The above equation is 
integrated numerically using Gaussian quadrature method as 
follows: 


m n 

j f(z)e dz = y H. f(a.) 
o i=l 1 1 


(4.24) 


where H^'s and a^’s are weight factors and zeros of nth order 
Laguerre polynomial. Typical values of ionization rate coeffi- 
cients for Hel are given in Table 4.2. 

For' collisional excitation coefficients k,(p r q),' the 
f oil© wing cases arise. fsr Hel system. 
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TABLE 4.2 

Values of He I ionization rate coefficients K(p,i) i=l in cm^sec 

at T g = 11.2 eV and 5.9 eV 



~ " ■ ” ‘ ‘ 

-- ' - - - 

- 

> - - - - -- ■ 

' - - - - 

Level* 

No. 

11.2 eV 

5.9 eV 

Level 

Mo. 

11.2 eV 

5.9 eV 

1 

_Q 

1.53. y 

1.34- 10 ** 

17 

4. 08" 6 

4.05" 6 

2 

-7 

2.12 

1.2 0~ 7 

18 

4.12" 6 

4. 09" 6 

3 

3.09 -7 

1 , 92 -7 

19 

6.42" 6 

6.69" 6 

4 

-7 

3.71 ' 

_7 

2.41 

20 

7 . 03" 6 

7 . 3 9" 6 

5 

-7 

2.27 

2. 87“ 7 

21 

7. 39" 6 

7 . 80" 6 

6 

1.22 -6 

1.0l“ 6 

22 

7. 5 8" 6 

8.03" 6 

7 

1. 47" 6 

1.26“ 6 

23 

7.59" 6 

8.04" 6 

8 

1.62“ 6 

1.41* 6 

24 

7.59" 6 

8. 04" 6 

9 

1 .7l“ 6 

1 . 5l“ 6 

25 

7.64" 6 

8.10“ 6 

10 

1.72“ 6 

1.51" 6 

26 

1.22" 5 

1.35" 5 

11 

1.74“ 6 

1.53" 6 

27 

1.82" 5 

2.07" 5 

12 

3 . 26 -6 

3 . 14“ 6 

28 

2.54" 5 

2.95" 5 

13 

3 . 68~ 6 

3 . 60~ 6 

29 

3.38" 5 

4.00" 5 

14 

3.89“ 6 

3.83“ 6 

30 

4. 36~ 5 

5.22" 5 

15 

4.08“ 6 

4. 04" 6 

31 

5.42" 5 

6.56~ 5 

16 

4. 08“' 6 

4. 04" 6 

32 

6.64" 5 

8.11" 5 

* 

Refer 

to Table 4. 

, 1 for the Hel 

states 




** —10 ~1 0 

Read 1.34 as 1.34 x 10 , etc. 
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*i 'll 

(i) Optically allowed transitions; (2 X P-1 X S, 3 P-2 S, etc.) 

The expression for excitation cross- sept ion q (E ) used by 

Pt- q e 

Drawin 50 is as follows; 

- e H —.2 

q (E ) = 4ir a 2 5 1 ’~ ■■ f(p,q) g(x) (4.25) 

qo, q e o ; E ^ 

where a is the Bohr radius, = 13.59 eV, E - E -E , 
o l pg q P 

f(p,q) is the absorption oscillator strength, E g the electron 

energy, and g(x) = « In (1.25 3 x), (here x is' •' 

x P,q 

equal to E /E )• a and 3 are constants close to unity. Using 

e P/ q 

47 

f (p,q) values from the Tables of Wiese and others and 
integrating the cross-section expression over Maxwellian 
electron energy distribution, excitation rate coefficients 
K(p,q) are obtained. 


(ii) Optically forbidden transitions without- change in •nrulti.-#. 
plicity: ■ (1 1 S-2 1 S, 2' 1 P-2 1 D, etc. ) 

50 . 

The excitation function given by Drawin is 


(E ) = 
q e 


4it a 




i (1 - A) 

q x x 


(4*26) 


whe ire x is defined as in (i) and Q is a constant for* a given 

P> q 
50 

transition and tabulated by Drawm. 


(iii) Optically forbidden transitions with change in multi- 
plicity: 

(a) Drawin 50 has given the following expressions for cross- 
section for transitions involving the ground state 
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(1 1 S-2 3 S / 1 1 S-2 3 P etc.) 


P/ q e 


(E ) = 4tr - 2 - x _1 


°P, 


q 5 ' 
^ x 


(4*27 ) 


and (b) those not involving the ground state (2 3 S-2 1 S / 2 1 S-3 3 P 
etc. ) 


q (E ) = 4ir a^ 
P/ q e o 


V 


x 


-1 


( 4 . 28 ) 


Qp^g values as given by Drawin are used. 

(iv) Excitation cross-sections which were not given by Drawin 

tor the case iii (b) were obtained from the compilations o£ 

51 

Moiseiwitsch. - As in (i), the rate coefficients are obtained 
by usual integration over Maxwellian energy distribution. 


Typical values of K(p,q) calculated according to the cross- 
sections above .are given in Table 4.3. 


The radiative recombination coefficient 3(p) valpes for 

Hel are obtained by appropriate integration of the recombina- 

. . 52 

tion cross-section given by Kramers : 


q 


(v) 


128 TT 4 Z 4 e 10 
3 ^ c 3 h 4 m v 2 n 3 


g(v) 


(4.29) 


e 

2 

where v is defined by hv = hv + m v 7 vis the frequency of 

o e 

radiation, v Q the threshold frequency and g(v) free-bound 
Gaunt factor, and Z is the nuclear charge, equal to 2. 

Typical values of B(p) obtained from eq. (4.29) are given in 


Table 4.4. 
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TABLE 4.3 

Typical values of He I collisions!- excitation coefficients K(p>.q:) ' 
and collisional deexcitation coefficients K(q,p) in cm! sec“l at 

T e = H.2 eV and 5.9 eV 



*r ^ v- ^ ^ 

- ■ ' * ~ ' - - 

- ' 

" * - ' -■ - ^ 

- ' 

p q 

11.2 eV 

5.9 eV 

p • q - -11-.-2- e-V - 

5.9 eV- 

l 1 S(l)-2 1 P(5) 

1.54- 10 * 

5 . 91 -11 

2 3 S(2)-3 3 P(8) 

_g 

5.30 * 

7 . 42” 9 


3.45“ 10 ** 

2.37“ 10 


2. 35~ 9 

-9 

2.07 * 

-3 1 P(ll) 

3.17" 11 

3.33" 12 

-4 3 P{14) 

1.29“ 9 

8. 11~ 10 


8. 42 -11 

5 . 73 -11 


6.11- 10 

5.24 _1 ° 

-4 1 P(18) 

1.26' 11 

i.36" 12 

-5 3 P(2l) 

5.42“ 10 

3 . 29 -1 ° 


3.82 -11 

2.59 -11 


2.64" 10 

2.24" 10 

2 1 S(3)-3 1 P(ll) 

2. 46 -8 

1.84“ 8 

2 3 S(2)-2 1 S(3) 

3.08" 8 

3.98“ 8 


1.03“ 8 

-9 

9.34 y 


9.93" 8 

1.37" 7 

-4 1 P(18) 

-9 

4.81 

-9 

3.30 * 

-3 1 S(7 ) 

6.43" 9 - 

6.9l“ 9 


2. 12~ 9 

-9 

1.88 y 


2.55“ 8 

3 . 52~ 8 

-5 1 P(25) 

1.59“ 9 

1.05" 9 

-4 1 S (13 ) 

5.01" 9 

5.07 -9 


7.20“ 10 

6.28 -1 ° 


2.13 -8 

2. 93~ 8 

l 1 S(l)-2 3 S(2) 

2.66 ^ 

7.14 -12 

2 3 S(2)-3 1 D(10) 

5. 67~ 9 

-9 

8.72 


-11 
5.27 J -- L 

6.98- 11 


4.56~ 9 

-9 

9.12 * 

-2 3 P(4) 

-11 

1.82 

4. 40” 12 

-4 1 D(16) 

-9 

2.03 

-9 

3.00 y 


1.33 -11 

1.74" 11 


_9 

1.73 y 

_9 

3.51 

-3 3 S(6) 

1.26 -11 

-12 

2.59 

-5 1 D(23) 

8. 76” 9 

_9 

1.23 


3. 23 -11 

. .,-11 
4. 16 


7. 68 _1 ° 

1.31 * 


n=6 (26) 

n=7 ( 27 ) 1.2l“ 5 

1 . 40“ 5 


-6 

-6 


4.47 

5.22 

n=ll(3l) 

n=12(32)l.70 -4 

-4 

2.09 


1.43“ 4 

•• ■1.76"’ 4 


* Upper values denote excitation coefficients K(p, q) 

**Lower values denote deexcitation coefficients K(q, p) 

The members in the parentheses represent the level mmber. 
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TABLE 4.4 

Typical values of Hel radiative recombination coefficients 8(p) 



in cm^ sec 

: “ 1 at T =11.2 
e 

eV and 

5.9 

eV 


evel No. 

11.2 eV 

5.9 eV 

Level 

No. 

11.2 eV 

5.9 eV 

Cp) 



'Cp)' 




l 

1.07" 12 

-1 2 

1.48 X 

17 


6.21“ 15 

8.56 -15 

2 

1.03" 13 

-13 

1.43 

18 


6.16- 15 

-15 
8.48 13 

3 

7.48" 14 

1.03" 13 

19 


3.7 8 -15 

-15 
5.21 ■ L3 

4 

6.38" 14 

8.72 -14 

20 


3.42 -15 

-15 
4.71 X3 

5 

5.64” 14 

7.74 -14 

21 


3.23 -15 

4. 45~ 15 

6 

2.14 -14 

2. 94 _1 4 

22 


3.14- 15 

4.32" 15 

7 

1.78 

2. 45~ 14 

23 


3.13- 15 

4.32 -15 

8 

1.62 -14 

2.23 -14 

24 


-15 
3.13 13 

4.32 -15 

9 

1.53 -14 

2.10 -14 

25 


3.11- 15 

4.28 -15 

10 

1.53- 14 

2 . 10~ 14 

26 


1.81" 15 

2.40" 15 

11 

1. 51~ 14 

2 . 07 ~ 1 4 

27 


1.13 -15 

1.55- 15 

12 

7.90" 15 

1.09 -14 

28 


7.53 -16 

1.04- 15 

13 

6.84- 15 

9.59" 1 ' 5 

29 


5.29" 16 

7.28 -16 

14 

6.54- 15 

9.02” 15 

30 


3.84 -16 

5.29^ 1$ 

15 

-15 

6.22 

8.57~ 15 

31 


2.91" 16 

4.01" 16 

16 

-15 

6.22 XD 

-15 

8.57 x 

32 


2.24 -16 

3. 09~ 16 
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8(i)/ the radiative recombination coefficient values 
for hydrogenic ion Hell are obtained from Seaton's 33 formula- 
tion 

6(i) = 5.197 x 10~ 14 Z x. 3//2 S’ (D) (4.30) 

i n 

where D = 157890 Z 2 /T and x. = D/i 2 . s (D) values are 

6 z n 

obtained by interpolation using Seaton's table. Typical 
values of f}(i) are given in Table 4.5. 

Rate coefficients for collisional excitation and ioniza- 
tion for Hell are obtained by using expressions used by Bates, 
Kingston and MCWhirter. 4 The Hell excitation coefficient is 


given by 

K(i,j) = - 4 * 7 ' 5 ' * 10 - 5 1 9 j j -v' 1 - exp (-E. AT J (4.31) 

Z j -i T^ 13 

where f^ is the absorption oscillator strength and E_^ = 

Ej - E^. Oscillator strength values are obtained from the 


table of Green, Rush and Chandler. 


48 


Hell ionization coefficients are given by 


vf . v 1.4 x 10 5 .2 , E ic v 

K(i,c) = - Y if ' — 1 ex P ( ) 

Z T| v e 


(4.32) 


where E^ c is the ionization potential of ith state. Typical 
values of K(i,j) and K(i,c) are given in Table 4.6. 

The collisional deexcitation coefficients K(q,p), K(j,i) 
and three-body recombination coefficient IC(i,p) and K(c,i) are 
obtained from eqs. (4. 12) -(4. 15) . 
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TABLE 4.5 

Typical values of Hell radiative recombination coefficient 


&(i) cm sec - 

at T = 11.2 eV 
e 

and 5.9 eV 

i 

11.2 eV 

5.9 eV 

1 

3-.29 -13 

6.94 -13 

2 

1.50 -13 

2.29 -13 

3 

7. 44 -14 

1.25 -13 

4 

4.12 -14 

7.69" 14 

5 

2.80 -14 

5.03 -14 

6 

1.92 -14 

3.63 -14 

7 

-14 

1.41 

3.53 -14 

8 

-14 

1.03 

2.09" 14 

9 

8.01 -15 

1.90" 14 

10 

6.26 -15 

1.22" 14 

11 

4.93 -15 

l.l -14 

12 

4.10 -15 

8.2 -15 

13 

3.33 -15 

7.4 -15 



77 


TABLE 4.6 


Values of excitation i<(i,j) 
of Hell at T 

e 


and ionization K(i, c ) coefficients 
11.2 eV and 5.9 eV 



11.2 eV 

u * . *ir* 

5. 9 eV 

K(i, j 

) cm 3 sec -3 


K ( 1 , 2 ) 

4.95 -10 

2.57 -11 

K ( 1 # 3 ) 

4.04 -11 

1.15 -12 

I<(1,10) • 

4.52 -13 

l/) 

rH 

1 

00 

CM 

• 

CO 

K(2,3) 

7.85“ 8 

5.90" 8 

K(2,4) 

8. 57 -9 

5.21 -9 

K(2 / 10) 

1.67" 10 

8. 10~ 11 

K(8, 9) 

1 . 79~ 5 

2.40 -5 

K(9,10) 

2. 79 -5 

3. 80~ 5 

K(i,c) 

cm 3 sec - ^ 


K(l,c) 

1.50" 10 

2.64“ 12 

K(2 / c) 

2.30 -8 

1.06“ 8 

K(5,c) 

3.98“ 7 

4.62“ 7 

K(10 /C ) 

1 . 85 -6 

2.44" 6 

K(l3,c) 

3. 18~ 6 

4 . 28 -6 
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Ionization- excitation coefficients K(p,i) for transi- 
tions Hel (p) -*■ Hell (i) are calculated using the following 


empirical expression . for excitation cross-section 


r—jjH - -*2 

fr (x) - 2 it a 2 ’ 1 ' ^(*) 

~ Z * a o ' E . r 

pi * 


(4.33) 


where g(x) = In (1.25 x), x = E /E . , and E .=E.-E . (4.34) 

x 2 e' pi' pi i p 


Ionization— excitation cross-section value q^ calculated from 

eq. (4.33) matches closely the cross-section presented for 

the case He(l) + e -*-He" r (4) +e+e by Drawin. It should be 

pointed out here that cross-sections for simultaneous 

Tohization-excita tion of Hell states from Hel ground state 

54 

was studied by Dalgarno and McDowell and their cross-sections 
are more than one order of magnitude lower than those predicted 
by eq. (4.33). Detailed analysis of the results given later, 
however, show that excitation function used here yield reasonable 
matching of theoretical and experimental results. Typical 
values of K(p,i) obtained from eq. (4.33) are presented in 
Table 4.7. 


Rate coefficients of excitation, deexcitation, ioniza- 
tion for processes involving collisions between the excited 
states of Hel and the ground state Hel(l) and the reverse 
processes, viz., KN(p,q), KN(q,p), KN($3,I), KN(l,p), are 
calculated using expressions given by Drawin, Klan and Ringler. ^ 
Typical values of these coefficients are given in Table 4.8. 



TABLE 4.7 


Values of ionization-excitation coefficients K(p,i) 

He (p ) -^Hell(i) cm 3 sec - '*' at T =11.2 eV and 5.9 eV 

e 


w w ' - s ' ' 

** - w * 

" - " 


(p/ i') 

11.2 eV 


5.9 eV 

1 1 S-1 

7.73" 10 


7. 2 -11 

l 1 S-2 

2.56 -12 


8.lO -15 

l 1 S-3 

-12 

1.03 


-15 

1.77 

1 1 S-12 

5 . 3 0” 1 3 


5.7 4 -1 ° 

3 

2° S-l 

1 . oo -7 


5 . 67 -8 

2 3 S-2 

3.26 -11 


5.18 -13 

2 3 S-3 

1.20 -11 


1.02 -13 

2 3 S-l 2 

-12 

5.80 


3.14" 13 

2 1 S-1 

_7 

1.47 


9.07 -8 

2 1 S-2 

3.63 -11 


6.16 -13 

2 1 S-3 

1.33" 11 


1.21 -13 

2 1 S-12 

6.4! -12 


3.70 -14 

n=10-l 

2. 54 -5 


2.94 -5 

n=10-2 

6.19 -11 


1.44 -12 

n=10-3 

2. 20 -11 


2.75~ 13 

n=10-i2 

l.OS -11 


8. 24“ 14 
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TABLE 4.8 

Values of neutral-neutral collision coefficients KN(p,q), KN(q,p), 

KN (p’ r i ) and ' KN ( %, p ) 

T g = 0.09 eV and T g = 11.2 eV, n(l) = 2x1 0 1 4 cm -3 


p 


q 

KN(p,q) cm 3 

— 1 3 

sec - KN(q / p) cm sec 

26 

(n=6) 

27 

(n=7 ) 1.38 -9 

1 . 20 -8 

26 

(n=6) 

28 

(n=8) 1.79" 10 

_9 

4.39 * 

26 

(n=6) 

32 

(n=12) 1.24 -11 

-9 

2.68 y 

30 

(n=10) 

31 

(n=ll) 5.15 -8 

8.14 -8 

30 

( n=l 0) 

32 

(n=12) 7.93 -8 

1.16 -7 

- 

P 

_ - — . 

KN(p,i) cm sec - 

KN (ijfp) cm 6 sec -1 




i=l 

' i=l 


25 

“ IT- 1 • • ' 

(5^) 

4.12 -14 

1.07 -29 


26 

(n=6) 

1.85 -11 

1.30 -24 


27 

(n=7 ) 

2. 38 -1 ° 

3. 48" 23 


28 

(n=8 ) 

8. 49 -1 ° 

1.36 -22 


32 

(n=12 ) 

1.73 -8 

4.57 -21 
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Optical Escape Factors : 

3 

Holstein has given the following simple expression for 
calculation of escape factors: 


A pq 


1.J5 

T(p, q) ( it In t )' 5 

pq 


(4.35) 


where x(p,q) is the optical depth. 

55 

Drawin and Emard have considered optical escape fac- 
tors in detail and included effects due to Stark broadening, 
inelastic and superelastic collisions suitable to collisional 

radiative model, and derived expressions for A and t 

P,q p, q 

1 1 

in terms of 2 P-1 S transition for Hel. But in the present 
study, electron density being low (10 10 -10 13 ) for Stark broaden- 
ing to become important, only Doppler broadening is taken into 
account and accordingly the following expression for optical 
depth is applicable: 


where 


and 


P/q 


_ T £ pq * 2 P.q 5/ Xd Pq 
12 f 1 2 (AX D^12 


(AA-) 

D pq 


12 


X 2kT . /n 
PQ , _ g-jl/2 

c V M 


m c 
e 


£ » nU) L & 12 


(4.36) 

(4.37) 

(4.38) 


Here (1,2) indicates the first resonance transition (l 1 S-2 1 P), 


f is. the absorption oscillator strength, x-' is the wave- 

pq pq 

length of the line radiation, T is the gas tanperature, M 

g 

and m e are the atomic and electronic masses, and L is the 
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average length through which the light travels. 

A lq values calculated from r (from eg. (4.35)) are 

used in the calculations wherein A(q, l) are replaced by 

A-, A ( q, 1 ) . 

1/q 

Plasma Parameters Used in CR Model Calculations: 

From probe data, axial variation of T g was found to be 

within 10% at all magnetic field values. Further, at high and 

low field values, the ion saturation current was found to be 

practically invariant of the axial position. At median 

magnetic field strengths the variation was found to be ^15% 

with respect to the value at the centre of the colunn — with 

high values at the electron gun end and low values at the pump 

end. On. the basis of these observations it was concluded 

that, although some axial variation exists, reasonable values 

of average population densities can be obtained using the n , 

T values at the centre of the column, 
e 

n Determination from Current Ratio., Method • 

e ......... 

In Chapter 3 results were presented on measurements of 
emission intensity ratios at pairs of emission currents for a 
number of He I lines. On the basis of these ratios and colli- 
sional-radiative model calculations, it is possible to derive 
electron density variations in the plasma as has recently been 
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a Ppli ec ^ to hydrogen plasma. ^ Examination of the r and r, 

o 1 

coefficients given in Table 4.10 as well as those in Ref. 43 
v/ill show that r^ coefficients, in the temperature and elec- 
tron density range of the present experiments, has negligible 
contributions to the population densities of excited states 
in comparison with contributions from r^ . Neglecting thus 
the contribution of r o in eq. (4.18), the rdtio of r^ values 

at two electron densities n and n is given by 

e 2 e 1 


n 


n 


r. "(p) n 2 (pi 


n 


n e 

I 2 (p) 
n 


(4.39) 


(p) n 1 (p) 


(p) 


n 


n 


X ^2 

where I (p), I (p) refer to omission intensities of a 
given, transition at the two electron densities. The second 
equality in eq. (4.39) is valid because of cancellation of 
the common transition probability appearing both in the 
numerator and the denominator. The above equation can be 
used for deriving the plasma electron densities. 

Table 4.9 shows that the ratio of r^ values at two elec- 
tron densities n andn is. dependent only on the magnitude 

1 2 

of electron densities and is insensitive to changes in electron 

h e n e 

temperature. Thus a certain ratio of r 1 2 /r 1 1 is unique to a 

certain n /n ratio and the n (or n ) value. If a plot is 
e 2 e l e l e 2 

made of this r. -ratio versus the electron density n (or n ), 

e l e 2 

and it is assumed that the plasma electron density n:'s are proportional 
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TABLE 4.9 
n e n 

Values of the ratio of r, Vr,* 1 for tho Hel 3 3 P state at 
electron temper 3 .tu.r 0 s .. 


n 

C-'< 

W 1 

n Q ratio 
= n e/ n e 

2 e l 

11*4 

T 

e 

9.5 

(ev) 

8.0 

5.8 

3.17 

1 . 0 10 

3 

2.98 

2.99 

2.99 

2.99 

2.99 


4 

3.97 

3.97 

3.98 

3.98 

3,98 

5.0 10 

3 

2.93 

2.94 

2.94 

2.95 

2.95 


4 

3.87 

3.88 

3.89 

3.90 

3.98 

l.O 11 

3 

2.88 

2.89 

2.89 

2.91 

2.92 


4 

3.77 

3.78 

3.79 

3.82 

3.96 

5 . 0 1 1 

3 . 

2. 62 

2.61 

2. 62 

2.67 

2.90 


4 

3.30 

3.29 

3.30 

3.39 

3.39 

l.O 12 

3 

2.41 

2.40 

2.41 

2. 42 

2.42 


4 

2.96 

2.92 

2.95 

3.04 

3.09 

l.O 13 

3 

1.51 

1.49 

1 . 49 

1.49 

1.50 


4 . 

1.62 

1.60 

1.59 

1.60 

1.63 
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to the emission current i *s (i as defined in Ch. 2), the 

e e 

electron density at which the r 1 -ratios match with the experi- 

j -e 9 i e 

mental emission intensity ratios I ^ (p) /I 1 (p) is the 

electron density of the plasma at emission current i (or i 1 . 

e l e 2 

By using emission intensity ratios- at emission currents 3 amp 

to 1 amp and also 2 amp to 1/2 amp, plasma electron densities 

were derived. The r^ ratio plots for the states 3^D, 4^D, 5^D 
3 3 

and 3 P, 4 P are given in Pigs. 4.1, 4.2 and 4.3. In these fig- 
ures rj/r^ (or r^/rj) represents the ratio of r values at 4n 

(or 3n ) to that at n - . Figures 4.1, 4.2 -and 4*3 in coniunc- 
1 e l ' - 

tion with Figs. 3.12, 3.13 and 3.14, are used to derive electron 

densities. Electron densities, derived in this manner versus 

the field plot is given in Fig. 4.4. r^-Ratio plots of 

Figs. 4.1, 4.2 and 4.3 can be made with any collisional- 

radiative model results, but in this case, r^ values obtained 

from the collisional-- radiative model calculations in the present 

work have been used. 

It is necessary, for calculation of population densities, 

to have knowledge about the plasma parameters n g and T . 

However, the n values deduced from the ratio method as above, 

0 

and n from probe data are at variance. The same is true for 

T e obtained from probe method and that from ratio of singlet- 

triplet emission intensities, as shown in Fig. 3.19. Rather 

than making collisional-radiative model calculations for all 

such plaana parameter conditions as shown in Fig. 3.19 and 

Fig. 4.4, a reverse approach is adopted, wherein n , T 

© 0 

profiles are gyolvied from concurrent considerations of the 
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observed enhancement factors and the predictions from the 

model. Specifically at B=0, n =4x1 0 10 , T =11.2 eV and at 

© 0 

B = 700 / n o = 3.2 x 10 73 , T = 5.9 eV are used. The n , T 
- e e e 

profile (as a function of the magnetic field B) adopted in the 
intermediate region is shown in Fig. 4.15. Later, such n , 

T e profiles are compared with those obtained from experimental 
measurements of plasma parameters and reasonableness of the 
Hg/ T^ profiles used in the calculations are discussed. 


Calculation of Enhancement Factors and Comparison with 
Experimental Results: 


As pointed out earlier, 43 equations involving the popula- 
tions of Hel and Hell were solved simultaneously, using IBM 

7044/1401 computer, to give final r , r, , r" r , r* and rt coef- 

o i o 1 2 

ficients from which populations of various Hel and Hell states 

could be calculated from eqs. (4.18) and (4.19). Calculations 

were made for an electron temperature range of 3 eV-14 eV at 

every 0.5 eV interval, an electron density range of 10 9 -10 17 

with 5 values (1,3, 5,7, 9) for every decade, neutral densities 
1 a 14 

of 10 to 4 x lO , at gas temperatures T = 500, 1000, 1500 

and 2000°K, and IIe + (l')/He ++ ratios !, 10, 1000, 1G 4 , lO 5 . 

Calculations were made for three sets of Hel optical escape 

factors: (i) for all radiation A = 1, i.e., plasma optically 

pq 

thin, (ii) plasma partially optically thiclc for resonance 


radiation with 5 = °* 000 9, Ag ^g = °* °°6, 


18 


= 0.014, 


A l,25 ~ °*° 3 ' A lf 


26 _ 0.07, 2? - 0.15, A 1/ 28 ~ °* 35; 
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the rest Apg = 1 , (iii) plasma partially optically thick 

for resonance radiation with a _ = 0.0095, A- = 0.075, 

■i / ^ 1 / J. J. 

*1,18 = 0 - 15 ' 5,25 = °' 35 ' 5,26 = °- 7 3nd rest *pq = ** 

(Note that p, q refer to level numbers given in Table 4.1.) 

The case (ii) corresponds to L ^ 40 crrm and 

(iii) to L ^6 cm in-eqi (4.38). For Hell, Lyman radiation (i) 
completely optically thick for all lines, i.e., A = 0 , 

(ii) A ^2 = 0 . 001 , A^^ = 0 . 01 , A^ = 0.06, A ^ 5 = 0.1 and 

rest Aj^j = 1 corresponding to L ^6 cm. Calculations were 
made for all possible combinations of n , T , ii(l),.. T~, h + )U )/n ++ , 
and various optically thick conditions detailed above. 

Some typical values of r , r^ are given in Table 4.10. 

For Hel(p) population density calculations, the ground state 

Hel(l) density is taken as ^2 x 10^. It will be observed 

that the population of Hel levels given by eq. (4.18) is 

dominated by the second term which contains the population 

coefficient r. , since contribution from r values are 
i o 

negligible in the entire range of electron densities involved 

in the experiment. r^ values for 1 set of conditions involving 

the five series of transitions are presented here, in Figs. 

4. 5-4.9. The conditions used are n (l ) =2xl0 14 , T '= lO0O°K, ; 

9 

T g = 5.9 and 11.2 eV, n + (l)/n ++ = 10^, and plasma partially 

optically thick Aj 5 = 0.0095, A^ ^ = 0.075, A-^ 18 = 0.15, 

A. oc = 0.35, A. = 0.7 and all other A =1. 

1,25 . 1,27 p,q 
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TABLE 4.10 


Values of r (p) 

O 

/ r- (p) 

coefficients and contributions 

from r (p) 
o 

and r^(p) to populatio 

n densities of some 

Hel excited 

states 

Partially optically thick case C A pq' 

s as in colunn 3 of T. 4.1l) 

J 

n n 

e 


T e = 11.2 eY, n( 1 ) = 

= 2.0 14 cm” 5 





r. (p7 n F (p") . 

n(l7 

»e P 

r o (p) 

^(p) r Q . n E (p ) 

rLXi) ' 

— • • n (p) =A+B 

^ -»■ 1 %rjti - k 

■>. ~ v - - r 

~_A 

SZj 

. . a. . 

. . , , Jem" 3 } ^ 

10 10 6 (3 3 S) 

5.5 

2.34“ 8 3 . 64 -3 

1.82 6 

1.82 6 

9 (3 3 D) 

3.2 

1.80" 9 5.6" 2 

6.19 S 

6.79 5 

' 11 O 1 ?) 

3.3 

2.10" 8 l.l -2 

1.58 6 

1.58 6 

15(4 3 D) 

3.8 

-9 -2 

1.04 y 6.3 ^ 

3.70 5 

3.70 5 

22 (5 3 d) 

2.3 

9.89" 10 2. 2~ 2 

3.40 5 

3 . 40 5 

10 14 6 

1.7 

1.78” 5 6.18 5 

g 

1.38 y 

1.38 9 

9 

1.1 

2. 34~ 6 1.93 6 

8. 82 8 

8.83 8 

11 

1.7 

2.28” 5 5.98 5 

1.72 9 

1.72 9 

15 

0.97 

6. 38 -7 1.6 6 

2. 27 8 

2. 28 8 

22 

• 0.98 

1.40“ 7 1 . 57 6 

4.83 7 

4. 98 7 



T e = 5.9 eV, n(l ) = 

2x1 0 1 4 cm” 3 


M 

O 

o 

ON 

0.36 

3. 29~ 8 3. 99 -3 

4.12 5 

4.12 5 

9 

2.00 

2.78~ 9 1.04 -1 

1.63 5 

1.63 5 

11 

0.14 

2. 48 -8 1. 42” 3 

2.91 5 

2.91 5 

15 

1.24 

-9 -2 

1.63 y 5.70 

8. 56 4 

8. 56 4 

22 

1.47 

-9 —2 

1.52 6.48 

7.57 4 

7.57 4 

10^ 4 6 

1.85 

2.42~ 5 2.03 6 

3.00 8 

3.02 8 

9 

1.14 

3 . 89~ 6 5. 93 6 

2.29 8 

2.34 8 

11 

1.88 

2. 25 -5 1.95 6 

2. 64 8 

2.65 8 

15 

1.00 

8. 9 1 -6 4. 65 8 

4. 68 7 

5.13 7 

22 

0.99 

1.7 8” 7 4.34 6 

8.86 6 

1.3.0 7 


...IV iA- *» ■- 









97 


The emission enhancement factor of a line, i.e. , the 
ratio of emission intensity at field B to that in absence of 
field is given, from a consideration of eq. (4.18) and 
neglecting r Q contribution, by the following equation: 


I B ( _ p) 

1'Cp) 


n B 


fP-, 

I Pi 0 

II 


and iP represent the 


exp (-E (- 1 ' - •-*-)) (4.40) 

P kT kT° 
e e 


B and that in the absence of field and E are as given in 

3? 

Table 4.1. Whereas in general such enhancement factors can 


be written for any two fields B^ and B 2 , we present here the 
enhancement factors wherein one of the intensities is at zero 


field. 


In order to evaluate such enhancement factors predicted 
by the model, it is necessary therefore to consider two r^ * 
plots, one for T at B=0 and the other for T at B=7 0O gauss. 
According to the discussion presented in the previous section 
we use T e (s=0) = 11.2 eV and T e (B=700 gauss) = 5.9 eV. This’ 
means that r^ plots at 11.2 eV and 5.9 eV are needed to 
evaluate the enhancement factors. However, if theoretical 
estimates of enhancement factors for the entire range of 
magnetic field variation is required, it is necessary to use 
the electron temperature profile as the magnetic field is 
varied and then use the appropriate series of r^ curves. Xp 
Table 4.11 the theoretical enhancement factors i b_7qc/ I B- 0 
axe presented for various conditions and then these are 
compared with the corresponding experimentally observed values. 
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TABLE 4,11 

Comparison of observed and calculated (Ig/I > 700 for 1 amp 

emission current 

Expt ” ' Theoretically calculated (ig/i^^Q 

Line Transition (ig/I ) 7Q0 Optically Partially* Partially** 

A fma-hlP T thin optically optically 

A =1 thick without thick with 

^ reduced reduced 


- -V 


_ 

1, , 

1 P-*- 1 D coeff. 1 P-^ 1 D coeff 
2 . _ 3 

5047 

4 1 S-*-2 1 P 

2.0 

4.7 

5.1 

5.0 

4437 

s^-s-^p 

2.0 

2.0 

2. 6 

2.0 

7065 

3 3 

3 S+2 P 

16.0 

13.0 

14.4 

13.8 

4713 

3 3 

4 S+2 P 

8.0 

9.3 

9.7 

9.2 

4121 

3 2 

3 &*-2 P 

2.0 

2.2 

1.9 

1.9 

5 875 

3 3 

3 D+2 P 

12.0 

12.7 

15.4 

13.2 

4471 

3 3 

4 D-*»2 P 

9.5 

10.3 

11.6 

9,8 

4026 

3 3 

5 D+2 P 

3.5 

4.4 

5.2 

3.0 

3 889 

3 3 

3 P-*-2 J S 

9.0 

8.5 

16.0 

10.0 

3187 

3 3 

4 P-»*2 S 

5,0 

4.7 

. 4.6 

4.5 

6678 

3 1 D -t-2 1 P 

14. 5 

10.8 

32.0 

12.2 

4922 

4 1 D+2 1 P 

9.0 

6.0 

19.5 

9.4 

43 88 

1 1 

5 D-*>2 P 

3.0 

3.0 

4.8 

2.9 

5015 

1 I 

3 P-*-2 S 

4.0 

17. *0 

6.6 

9.0 

3964 

4 1 P*‘2 1 S 

3.0 

13.8 

- 6.0 

6.0 ■- 

— 

- - 

- 

— ~ *■ ' 

• — 

- 

* V 

5 = 0, 0009, 

A l,ll = 

0.006, A, 

a t 

18 = 0.014, 

A l, 25 = °*° 3 

V 

26= ^ 

At, 27 = °* 

15/ A l, 28 

= 0.35, Rest 

A =1. 

pq 

** A 

1 

c = 0.0095 

s 

M 

1 — 1 

If 

= 0.075, A 1 

,18 = °- 15 ' 

A l,25 = °* 35 ' 

A 1 

,26 = °* 7 ' 

Rest A 

pq 

= 1. 
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For these calculations n (B=0) = 4 x 10 10 , n (B=700) = 

6 0 

3.2 x 1 0 1 2 i t = 1000 K, n + (l)/n ++ = lO 4 . 

9 


In Table 4.11, the first colunn of theoretical values is 
for the optically thin condition. It will be noted that the 
theoretical values of the ratio is within 9% except for the 
transitions 4 1 S-»-2 1 P, 3 1 r>2 1 P, 4 1 D-*2 1 P, 3 1 P+2 1 S and 4 1 B+2 1 S / 


out of which the last two lines show very large variations. 
Criteria for application of optically thick condition have 
been discussed in the previous section, and such condition is 
applicable to the results of the present experiments. Results 
of calculations using partially optically thick condition with 


Vs = °- 0009 ' "l,ll = 0.006, A 1j18 = 0.014, A 1i25 = 0.03, 

"1,26 = °* °7 , 27 = °-15, % 28 = and - t ^ le rest tran- 

sitions a = 1 are shown in the second colunn of theoretically 
P9 

predicted results. 


The immediate effect of application of optically thick 
condition is to increase the populations of 3 P and 4 P at low 
electron densities without significantly affecting the popula- 
tions at high electron densities, and this leads to a lowering 

of emission intensity enhancement factors. But this condition 

1 1 

simultaneously increases the populations of the 3 D and 4 D 
states for moderate and high densities and hence very large 
enhancement factors result. In this case the average varia- 
tion of the theoretical Ig/I Q values with respect to the 

1 1 

experimental ratios is within 15% except for the lines 4 S+2 P, 
and 4 1 D+2 1 P. 
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While exploring possible conditions which might further 
bridge the gap between the experimental and the theoretical 
ratios, one notes that use of optically thin condition and 
higher values of collisional excitation coefficients from 1 1 S 
to n P would also lead to lower values of emission enhancement 
factors for the lines 3^P-*2^S and 4^P-*-2^S. However, not only 
does this not take care of the requirement of partially 
optically thick condition as previous discussion indicates, 
but the rate coefficients of the resonance transitions 1 7 S— n^£ 
are quite well studied, and the values seem to be reasonably 
accepted. On the other hand, if partially optically thick 
condition is used as shown in the column 2 of the theoretical 
results, additional considerations are necessary which will 
at least lead to lower values of the enhancement factors for 
the lines 3^D-*-2^P arc’ 4^D-*2~P. In this regard, two changes 
can be incorporated, primarily with the idea of finding condi- 
tions which makes the theoretical and experimental ratios 
match closer, (i), which involves making the plasma less 
optically thick compared to that used for calculations of 
column 2, and (ii), which uses reduced values of the excitation 
rate coefficients of ^P-^D transitions. Specifically, the 
plasma is considered partially optically thick with A =0. OQ95, 

^1,11 = °*° 75 ' *1,18 = °* 15 ' A i,25 = °* 35 ' a i, 26 = °* 5 and 

the rest A pg = 1, and the K(p, q) used for 1 P- 1 D transitions 

are 1/2 of the values that can be calculated using eq. (4.25). 
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The results of such calculations are shown in column 3 of 
theoretical results in Table 4.11. Among the three sets of 
theoretical results presented, the last set of calculations 
given in column 3 provides the best match with the experi- 
mental enhancement factors giving an average deviation or 
1 %, except the lines, 3 1 P-»-2 1 S and 4 1 S-»-2 1 P. 

While Table 4.11 shows the emission enhancement fac- 
tors l(B=700)/l(B=0), emission enhancement profiles for the 
entire range of magnetic field variation, calculated by 
using the conditions of column 3 of Table 4.11, and n e , T^ 
profiles of Pig. 4. IB are presented in Figs. 4.10-4.12 for 
the transitions 3^D*»2^P, 4^D-*-2^P and 5 3 D-»-2^P, 3 P + 2 S and 
4 3 P-K2 3 S, 3 3 D-K2 3 P, 4 3 D-»-2 3 P and 4 3 D-*2 3 P. 

It is to be mentioned here that the fit of the theoretical 

and experimental enhancement profiles is not significantly 

inferior from Figs.- 4.10-4.12 if partially optically thick 

3 3 

conditions of column 2 are used, except for the lines 3 P-^2 S, 
3" l 'D-^2 1 P, 4 1 I>*-2' ! 'P. This, however, shows up also in Table 4.11. 

All the experimental results presented in Table 4.11 
are for an emission current of 1 amp. Enhancement factors 
were also measured at emission currents of 1/2 amp, 2 amp 
and 3 amp. The experimental emission enhancement factors 
have already been presented in Table 3.1. It has been shewn 
that the electron densities are proportional to emission 
currents from data on electron collection current in the plasma 
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Fig - A • 10- Comparison of experimental and theoretical 
emission enhancement profiles 



o jgp 
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Fig. 4.12 Comparison of experimental and theoretical 
emission enhancement profiles 
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column, and also from probe data. The ion saturation .current 
was in the ratio of %: 1: 2: 3 for emission currents of % amp: 

1 amp: 2amp: 3amp. Theoretical emission enhancement factors 
for other emission currents therefore can be calculated 
using electron densities in proportion to that used for 
calculations at 1 amp of emission current. Electron tempera- 
tures are also necessary for calculation of populations. In 

longitudinal magnetic field electron temperature may be 
39 

considered to be independent of electron density. Experi- 
mental relative intensity measurements using He lines 4^S/4^S 
and 4 D/4 D to obtain T , as shown in Fig. 3.17 and 3.18, 
bear this out. Though the absolute values obtained are 
somewhat different, it is clear that at a particular B, T e ' 
does not have n g dependence in the experimental n g range. 

One, therefore, may use the same T^ range for 2, 3 amp of 
emission currents as has been used for 1 amp of emission 
current. Results of such theoretical calculations of enhance- 
ment factors i(B=7 00)/i(b= 0) are shown in Table 4.12. Some 
enhancement profiles are shown in Fig;' 4.13 which ate 
evaluated in the same manner as enhancement factors at 1 amp 
of emission current. It must be mentioned that, here also, 
when the results are compared to -the profiles corresponding to 

the conditions of go limn 2 of Table 4.11, the major variation 

3 - 3 . i ■ -i 

from experimental curves arises for the lines 3 P-2 Si 3 D~2 P 
and' 4 1 D-2 1 P. 
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TABLE 4.12 

Comparison of experimentally observed and calculated (Ig/l ^ ) 700 


for 0.5 

amp, 2 

amp and 

3 amp 

emission 

currents 


Emitting 

State 

Line 
% .. 

i e 

Cal ' 

0 . 5 amp 

Obs 

Cal 

amp 

Obs 

^ = 3 
Cal 

amp 

Obs 

3 3 S 

7065 

20 

32 

10.7 

10.5 

7.4 

8.0 

4 3 S 

4713 

14 

17 

6.2 

6 . 0 . 

3.9 

3.5 

5 3 S 

4121 

2.8 

3 

1.5 

1.7 

1.1 

1.2 

3 3 P 

3889 

14.0 

18.0 

7.3 

7.0 

5.2 

5.0 

4 3 P 

3187 

10.04 9.5 

3.6 

3.9 

2.2 

3.5 

3 1 D 

667 8 

16.2 

24 

10.5 

12.0 

9.5 

9.0 

4 1 D 

4922 

13.0 

14.0 

4.4 

4.8 . 

3.6 

3.0 

5 1 D 

4388 

5.2 

4.8 

2.5 

2.3 

2.1 

1.8 

3 3 D 

5875 

16.0 

17.5 

10.5 

9.5 

8.0 


4 3 D 

4471 

12.7 

12.0 

7.5 

6.0 

5.2 


5 3 D 

4026 

6.0 

4.5 

2.7 

2.5 

2.1 


4 1 S 

5047 

6.2 

2.8 

2.2 

1.5 

2.0 


5 1 S 

4437 

3.0 

2.6 

1.7 

1.5 

1.5 


1 

3 P 

5015 

10.8 

8.0 

7.5 

3.4 

5.7 


1 

4 P . 

3964 

8.0 

4*Q 

4.8 

2.3 

3.0 
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Fig *4.1 3 Comparison of experimental and theoretical 
emission enhancement profiles 
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Population of Hell levels also can be determined in 
a manner similar to that of Hel, i.e. using r*, r^ and r* 
coefficients (obtained by simultaneously solving the Hell 
equations with Hel ones) and assuming Hell(l) = n £ and Hel(l) 

^ 2 x 10 14 . Typical values of r*, r^ and r^ coefficients 
and some population values are shown in Table 4.13. Enhance- 
ment factors I (B=7 00)/I (B=0) are also shown in Table 4.13. 
Enhancement profile for i=4/ calculated by- using the T g and 
n e profiles of Fig. 4.15, in the range B=0 and B=700, is shown 
in Fig. 4.14. 

Now that the experimental and calculated enission 

profiles have been compared using the n , T profiles of 

Fig. 4. 15, it is proper to discuss these plasma parameter 

profiles in the perspective of experimental plasma diagnostics 

data. Electron densities measured using probes is found to 

be 7-8 times lower than those observed from current ratio 

method although the pattern of variation from low field to 

high field is the same. Although the probe data gives an 

order of magnitude check of n g , the probe n 0 is likely to be 

less reliable for 'density determinations' due to uncertainties 

of "-the .probe theory in the magnetic field. Studies 'of -Shaw 
57 

and Kino show that ion collection at the probe is normally 
lower in a magnetic field, corresponding to the plasma qf 
the same density in absence of the field, and point out that 
it is difficult to get complete ion saturation. That the’ -.adopted 
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TABLE 4.13 

r*(4), r*(4) and r*(4) of Hell (i=4) at 11.2 eV and 5.9 eV 


n 


r + (4) 

o 

r+(4) 

+ 

r 2 

(4) 

s 

11.2 

eV 5 . 9 eV 

11.2 eV 

5.9 eV 

11.2 eV 

5.9 eV 

io 10 

0.34 

0.20 

2.94" 8 

2. 06 -1 ° 

5 . 7 -1 ° 

6.5- 10 

10 12 

0.36 

0.23 

2.92" 6 

2.05 -8 

5.6- 8 

6. 2 -8 

10 14 

0.60 

0.40 

2. O -4 

' 1 . 8 -6 

4.0" 6 

4. 8~ 6 


T e eV 

n g Optical thickness 

n(l) 

n + (4) 

11.2 

4.0** 8 ' Optically thin 

14 

2.0 

9. 97 8 

5.9 

3.2 42 * Partially*optically 

14 

2.0 

4 

4.84^ 


thick 




Observed (I B /I 0 > 700 = 

C 1 ciliated 

^700 ~ 


6.3 

,4.84 

9.97 


at 1 amp 
4 

v = 4.85 


* A 12 = 0.001, A 13 = 0.01, A l4 = 0.06, A 15 = 0.1, & A = 1. 
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n e profile in Fig. 4.15 is located higher than the probe n g 
profile, is compatible with these observations. An inspection 
of Fig. 4.15 shows the - remarkable similarity between probe n e 
profile and the n Q profile used in calculations as well as 
that obtained from emission current ratio method, it may be 
noted that if the effective collection area is lower than the 
geometric area used in n 0 evaluation, significant improvement 

in n will result. 

© 


The electron temperature profile adopted for collisional- 

radiative model calculations is closer to T^ profile obtained 

from electric probes but considerably off from T g profile 

obtained through singlet-triplet ratio measurements. According 

3H 

to Bickerton and von Engel, electron temperatures measured 

by probes in the presence of magnetic field is almost at 
good as in absence of the field when the negative probe charac- 
teristics are taken for T estimation. They point out that 
any error in evaluation of probe T 0 is essentially due to 
assumptions inherent in the probe theory applicable to plasmas 
in absence of magnetic field. The presence of magnetic field 
does not alter the situation substantially. 

Electron temperature measurements from line intensity 

41 

ratio method has considerably more limitations. The ordi- 
narily used expression of singlet-triplet intensity ratio 
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(eq. 3.2)* is considered to be independent of electron density. 

However, if one computes the ratio of I 5047 /X 4713 various 

T' e and n Q values using a detailed collisional-radiative model/ 

one observes considerable dependence on n g at least in the 
11 13 

range 10 tolO , the n g range of the present experiments, 

though at higher and lower densities the ratio is essentially 

independent of n g . This is shown in Pig. 4.16. Hence at 

intermediate densities T obtained from this ratio method is 

e 

likely to be less reliable. 

It may be mentioned here that in the present work an 

attempt was made to obtain T g estimates from Boltzmann plot 

3 3 

for higher lines of n D-*2 P series but the intensities of the 

lines being very weak, measurements could not be made, it Has - 

been indicated however that this method gives too low T 
4*1 

values. 

C O 

It may be added here that Johnson and Hinnov J have 
pointed out that no acceptable line ratio method for electron 
temperature determination has been found in the density range 
of 10 10 to 10"^ cm - ^. They further comment that absolute 
intensity measurement of a line should be reasonably good 
for T measurement. But T obtained by them using a comparison of 
the absolute intensities - observed with those calculated from a, 

CR model for He plasma is reliable within only <*30%. 

Considering all the above .facts, it may be taken that 
the electron temperature obtained by using probes is probably 




cuiated emission inter 
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closer to the actual situation. With respect to the T g profile 
obtained from probe, the profile that gives best match with the 
experimental" enhancement factors has a maximum deviation of 20% 
at intermediate points end 5% at the end points. 

Effect of inclusion of more levels in the modei on population 
densities: 

Calculated population densities taking energy levels of 
He I up to principal quantum number n=10 and Hell levels up to 
n=13 were compared with those obtained by taking n up to 12 
keeping Hell levels n=13. In the n , T q range of the present 
investigation, maximum variation found in the calculated 
population densities for any excited state was less than 1%. 

It may be recalled that in Johnson 1 s^' model n=25 "was used. Blit 
the present model with n=12 matches with Johnson's experimental 
results when appropriate n Q , T , n(l) are- used, as does Johnson's 
n=25 model. From this, it seems that additional levels are 
quite unnecessary. More detailed comments on Johnson's experi- 
mental results will be presented later in this chapter. 

Effect of neutral-neutral collisions on the population densities: 

Variation in the population densities by including 
neutral-neutral collisions when compared with those in absence 
of such process was less than 2%. From an inspection of the 
values of KN(p,q), KN (p, i) from Table 4.8 and seme of the 
electron impact rate coefficients from Tables 4# 2 and -4.3, 
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such small variation seems to be reasonable because the electron 

impact collisions are dominant at T ranae of 12-5 eV. 

e 

Neutral He temperature T : 

g 

Neutral temperature estimate was made from Doppler widths 
of some of the He I lines. The temperature obtained was about 
300° K. Variation of this temperature in the magnetic field was 
not more than 20%. Hence an upper limit of T g equal to 100O°IC 
is used in the calculations for enhancement. However, calcula- 
tions in the T^ range 500-2000°IC does not make any noticeable 
difference in the results. 

Effect of variation of n g , ranges on Ig/l^Ccalc) : 

In order to find out whether any other range of n g and 

values would give better correlation with the experiments, 

calculations of Iq/I q were made chosing different sets of n e 

and T q values. Keeping the variation in T g same (11.2 to 5.9 

eV) and also (n"Vn°)„~„ = 80 as is taken in all calculations, 

0 © / 

l-n/I at 70C is calculated with n° = 10^°. The results show 
B o e 

that calculated I q /I q values are 40% higher than the observed 

values for first and second line of any series, and for higher 

o 11 

lines error exceeds 80%. With n =10 , calculated values 

e ... 

are found to be about 20% lower than the experimental 

values and here again the error is more than 100% for higher 

levels. Keeping n° = 4 x 10 i0 and T ranging from 12 to 3 V 

^ 0 



117 


(as observed in line ratio method ), it is found that a 10 3 
rold increase in the electron density is required to make the 
calculated Ig/I Q agree with the experimental results. Such 
large increase in.n is not supported by the present experiments. 

Absolute population densities of He + (4) : 

Although in this work absolute population densities of 
any level has not been measured, dependability of the calculated 
population densities can be checked by comparing the calculated 
values with those experimentally observed by other workers. It 
will be shown, later in this chapter, that the present model 
predicts absolute population densities accurately within a 
factor of 2 of some heliurt plasma experiments of other workers 
where absolute densities were experimentally measured. If this 
is also assumed to be valid in the present work, then from the 
experimental relative intensities of Hell 4-»-3 and Hel 4 3 Sw2 3 P • 
lines (0.09 at B=0 and 0.059 B=700 gauss), n(4 3 s) calculated from 
the model, and known transition probabilities of -the lln.es, one 
obtains for n + (4) = 1.8xl0 4 (at B=0) and n + (4) = 1.04 x 10 5 
(at B=700 gauss). If only heliurt ion model is considered such 
high population densities cannot be explained regardless of the 
degree (0-100%) of second ionization one assuttes. This seems to 
strongly indicate that Hel states must be regpcsnsible in a major 
way for populating the Hell excited levels. 
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Effect of Tt + (b)/n ++ 


variation on Hel and Hell populations : 


In the n g , T g range of the present study, Hel(p) states 
are populated mainly from Hel ground state and as a consequence, 
variation of n" f "/n ++ ratio has little effect on n(p) values-; 
since the n /n ratio would influence n(p) only when contri- 
butions from r (p) to n(p) become significant. In the calcula- 
tion of Hell coefficients r*(i), r*(i), and r*(i), from eqs. (4.16) 
and (4.17), n + /n ++ appears as a parameter in the term 
P (p) ,n e . K(p, i ) . ng (p)/n^"(i) and explicit values of the n + /n + " r 
ratios are used in solution of the simultaneous equations. 
Rewriting eq. (4.19), population densities of Hell(i) are 
given by 


,-16 


n + (i) = r*(i) n e .n ++ . ^ exp + *^(i) n(l) 


++ 

- n ' exp (ET - E )/kT 
n + (l) 1 1 e 


+ r^Ci) . i^ n^tlKexp (E^ - Eip/klg 


* 2 


(4.41) 


where E^ is the ionization potential of the Hell state i and 
is the ionization potential of the ground state of Hel to the 
first continuum. 

For a given n , T , n(l) and n + /n ++ values, n + (i) can 
0 0 

be calculated from eq. (4.4l). Whereas r + (i) and r*(i) are 

o z 

| 

functions of only n and T , rT(i) is a function of n , T and 

e e l e e 

n + /n ++ . At T >8.6 eV, n < 1(4 ^ and n(l) ^1<4^, major contri- 
0 . 0 

bution to n + (i) arises from the second term in eq. ^4.4l), while 
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- 19 4 - 

at T e <8.5 eV and n g > 10 major contribution to n (i) is from 

the third term. At very low electron temperatures and high 

densities (T <2 eV, n >10^), the first term in eq. (4.4l) 

0 0 

becomes important. 

Taking n g = n + +2n' + ,t it follows from eq. (4. 41 ), that if- the 
contribution to n^Ci) is from the first or the third term, varia- 
tion of n + /n ++ changes the n + (i) values to an extent determined 
by the variation in n ++ or n + (l). For T q >8.5 eV, it is found 
by calculation of n + (i) that its (viz., n + (4) etc.) variation is 
not significant when the parameter n /n is varied. Hot only 
does this show that Hel states are pr.iniarly responsible for 
populating the Hell(i) states, but also the dominant contribution 
originates from Hel ground state. This is because, a comparison 
between the metastable and the ground state ionization coeffi- 
cient K(p, i)*' s show that the value for the metastables is about 
1-2 orders of magnitude higher (Table 4.7) but the population 
densities are at least 3 orders of magnitude lower leading to 
relatively .insignificant contribution to -Hell (i) from the 
metastables. From higher states of Hel, the contributions are 
even less. 

* . . ■ 

Effect of optical thickness on the population densities : 

The extent of optical thickness for Hel radiation does 

. -f 

not significantly affect- He levels and conversely the optical 
thickness for He + radiation does not alter Hel populations much. 
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In general, if the plasma is optically thick, population densi- 
ties of all the excited levels increase, compared to the optical 
thin case. As a result of photo absorption the increase in 
population densities of the states directly, populated due to 
absorption (2-.P13 P ete.) is largest. For other levels 

which are collisionally connected to these states the increase 
is less. In the He I system, at lower electron densities, 

absorption of resonance radiation increases the population of 

l A 

n P levels to an extent there is decrease in^A(p, l) values. 

If the plasma is partially optically thick then, increase in 

population densities of the singlet states is higher than the 

corresponding states in the triplet series. 

The difficulties in ascertaining correct values of optical 
escape factors applicable to a particular plasma is well recog- 
nized. This arises mainly due to the uncertainties in L in 
eq. (4.38). In addition, even small .variations in neutral 
density significantly alter the escape factors. However, the 
fact that for optically thin conditions triplet level enhance- 
ments are explained quite well and n^P level enhancements are 
not, is a strong qualitative indication that the population of 
■^P levels involve some special processes. Reabsorption of 
resonance radiation, i.e., partially optically thick condition, 
seams to be the most reasonable process one can conceive of even 
though the exact values of optical escape factors may have 
uncertainties . 
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Collisional- radiative Models of other Workers Applied to the 
Present Experimental Results: 

Having presented a comparison of the experimentally 
observed Hei and Hell emission enhancements with the theoretical 
predictions of the model used in the present work, it should be 
be of interest to consider calculations of other workers on the 
helium system and check to what extent the present experimental 
observations can be explained with their theoretical results. 
First, calculations of Drawin and Emard^ are considered. 

Using Table 2 in Ref. 43, enhancement factors were calculated 
and are presented here in Table 4.14. The major difference 
between Drawin and Ehiard's calculations and the present one is 
that, in the present case, the shblevels of n = 3,4,5 HeX states 
are completely separated. In Table 4.14 the experimental values 
of emission enhancements and those calculated in the present 
work are also given. One observes that the present calculations 
provide much better correlations with the experimental enhance- 
ment factors of the sublevels than the calculations of Drawin 
and Emard using coalesced levels. It may be noted that popula- 
tion densities of sublevels obtained from those of the coalesced 
levels using appropriate Boltzmann factors would not give any 
better enhancement values than shown in Table 4.14. This is due 
to the fact that the difference in energy levels between the 
sublevels is very small even in the case of n=3 state. n=4 and 
n=5 states have still smaller energy differences between the 
sublevels and the situation would not improve. 
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TABLE 4, 14 


Comparison of emission enhancement factors (ig/I ^ ^ -jqq calculated :: 
from Drawin and Ernard*s tables (Ref. 43) with the present calcula- 
tions and experimentally observed values 


Emitting 

Triplet levels Calc. (Ref . 43 ) 


n=3 14.0 

(3 3 S+3 3 P+3 3 D) 


n=4 5.6 

(4 3 S+4 3 P+4 3 D+4 3 F.) 


n=5 2.25 

(all triplet levels) 

Emitting 
Singlet JLevels 

n=3 7.5 

(3 1 S+3 1 P+3 1 D) 

n=4 5.85 

(4 1 S+4 1 P+4 1 D+4- 1 F) 


n=5 

(all 


singlet levels) 


Experimental 
observation 
present work 

Calculated 
val ues 

present work 

3 

3 S 

16 

13.8 

3 3 P 

9 

10.0 

3 3 D 

12.7 

13.2 

4 3 S 

8.0 

9.2 

4 3 P 

5.0 

4.5 

4 3 D 

9.5 

9.8 

5 3 S 

2.0 

1.85 

5 3 D 

3.5 

3.0 

3 1 D 

14.5 

12.2 

s 1 -? 

4.0 

9 . 0 

4 1 D 

9.0 

9.4 

4 1 S 

2.0 

5.0 

4 1 ? 

3.0 

6.0 

5 1 D 

3.0 

2.9 

5*S 

2.0 

2.0 


2.9 
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59 

Turning to Park's work on coliisional excitation- 
deexcitation rate coefficient calculations, one notes that in 
his scheme of Hel energy levels, he has used separated levels 

3 

up to n=3 but uses combined singlet- triplet groups of 3 D 
and 3 1 D, 4 3 S and 4 1 S, 4 3 P and A 1 ?, 4 3 D and 4 1 D, 4 3 P and 4 X F, 

5 3 S and 5^S, 5 3 P and 5^P, 5 3 D and 5^D, a groining of the 
levels 5 3 F, 5^F, 5 3 G and 5^G, and fully coalesced sublevels 
for n=6 to n=20. Talcing this scheme of levels and the rate 
coefficients provided by him, and using other atomic para- 
meters as described earlier in this chapter, a computer 
programme was developed to calculate Hel level populations 
from which enhancement factors could be obtained. Hel 

V 

enhancement factors calculated from Park’s model as well as 
experimental values from the present work are shown in 
Table 4.15. It can be seen that the deviations are large. 
Part of this reason seems to be in his using coalesed triplet 
and singlet states; discrepancy arises because the triplet 
and singlet level populations change in markedly different 
manner with electron density and electron temperature. Also, 
in his model, l-jn' excitation for n >6 are neglected. 

Applications of the Present Model to Experimental ' Results 
of Other Workers : 

Though the absolute values of population densities 
obtained in present calculations seem reasonable, one of the 
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TABLE 4.15 

Comparison of ( Ig/ I Q ) 7 00 values obtained by using Park's (Ref. 59) 
energy level scheme and his rate coefficients with those obtained 

in present calculations 


Emitting 

State 

^ I B //I 0'*700 calcu 
lated using Park's 
scheme 

{ I /I ) * 

x B / 0 700 

calculated using 

present scheme 

^ I B //l 0^7OO 
experimental ly 
observed 

3 3 S 

11.73 

13.8 

16 , 

4 3 S 

1.87 

9.2 

8 

5 3 S 

1.0 

1.95 

2.0 

4 1 S 

1.87 

5.0 

2.0 

5 1 S 

1.0 

2.0 

2.0 

3 3 P 

8.0 

10.0 

9.0 

4 3 P 

13.39 

4.5 

5.0 

3 1 D 

3.7 

12.2 

14.5 

4 1 D 

10.2 

9.4 

9.0 

5 1 D 

7.2 

2.9 

3.0 

S 1 ? 

6.41 

9.0 

4.0 

4 1 P 

13.39 

6.0 

3.0 

3 3 D 

3.7 

13.2 

12.0 

4 3 D 

10.2 

9.8 

9.5 

5 3 D 

7.2 

3.0 

3.5 



* Values corresponding to Table 4.11, coltamn 3. 

** Values corresponding to Table 3.2, for i =1 amp. 
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limitations of the present work is that absolute population 

densities could not be obtained from the experiments, due to 

lack of a standard lamp in the laboratory and its not being 

available in the country. While such deficiencies can be 

made up later with the development of the laboratory facility, 

it is of interest to try out the model, which seems to explain 

reasonably satisfactorily emission enhancement factors, to 

predict absolute population densities in experiments of other 

1 3 

workers on heliun plasmas. For such a comparison Johnson's 

experiments in Princeton stellarator C is chosen. His experi- 

12 13 

mental condition involve n rancre 4.2x10 -3.3x10 , T 

e e 

range 14.4 eV - 4 eV, and neutral density 4x1 0^ ^-S. 2x10^ . 

For some of these conditions optically thin model can be used 

in a straightforward manner, but with increasing neutral 

density, application of partially optically thick condition 

is required. While calculations were carried out for the 

entire range of his experimental conditions, and matching is 

satisfactory, the results for some of the conditions are 

presented in Table 4.16 and Table 4.17. Table 4.16 presents 

cases where optically thin conditions can be used. In Table 

13 

4.17, at the highest neutral density of 8 . 2 x 10 , partially 

optically thick condition is used. Calculation of the optical 
thickness was made by taking 3.5 cm as the radius of the plasma 
cross-section (the maximum distance the light travels from the 
centre of the plasma) at 50 kilogauss and ecTti ( 4 . 35 ) applied 
to calculate the optical escape factors. Examination of 
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TABLE 4.16 

Comparison of Johnson's experimentally observed values of He I 
population densities (Ref. 13) with present calculations using 

optically thin condition 






— ' 

» - - 

- - - ' 



Hinnov 

' Present 

Hinnov 

Present 

Hinnov Present 



expt. 

calc. 

expt. 

calc. 

expt. calc.— 


T eV 

10.8 

a. 

5 

7.2 


0 

n e 


4.2 12 

1. 

9 13 

1.49 13 

States 

n(l ) 


8.8 11 

6. 

3 12 

8.5 12 

4 1 S 


6. 6 5 

4.79 5 

2. 2 6 

2. 6 6 

1.48 6 2. O 6 

s^-s 


2.4 5 

1. 24 5 

, .5 

6. 4 

6.06 5 

3.3 5 4.3 5 

a 1 ? 


6.9 5 

5 

3 . 0 3 

7.4 6 

4.5 6 

4.5 6 3.1 6 

4 1 P 


3.3 5 

2.02 5 

1.96 6 

2.4 6 

1.17 7 1.6 6 

5 1 P 


1.24 5 

1.06 5 

5.0 5 

7. 65 5 

5 5 

4. 8 3 5.1 * 

3 1 D 


5. 1 5 

3.28 5 

4. 0 6 

3.57 6 

2.4 6 2.7 6 

4 1 D 


3 . 0 5 

1.64 5 

1.45 6 

1.23 6 

9. 1 5 9.2 5 

5 1 D 


9.5 4 

6. 0 4 

4.3 5 

2.78 3 

5 5 

2.9 2.2 

3 3 S 


2.7 6 

1.35 6 

1. 2 7 

1.6 7 

7 7 

1.06 1.4' 

4 3 S 


7 . 3 5 

5.6 5 

2. 6 6 

4. 0 6 

1.87 6 3 . 4 6 

5 3 S 


1. 86 5 

9.55 4 

5 . 8 5 

6. O 5 

4.2 5 5. I 5 

3 3 P 


2. I 6 

2.42 6 

9 . 7 6 

2 . 0 7 

7 . 9 6 1 . 6 7 

4 3 P 


4. 9 5 

4.99 5 

1.82 6 

2.6 

1.40° 2.2° 

5 3 P 


1.49 5 

4 

9.0‘ 

4.7 5 

. .5 

4.4 

5 5 

3. 6 3 3. 6 3 

3 3 D 


1.26 6 

7.7 5 

7.8 6 

9.4 6 

5.0^ 7.7^ 

4 3 D 


5 

3.3 

5 

3.11 3 

1.53 6 

2.4 6 

1.13 6 2.0 6 

5 3 D 


5 

l.O 3 

1.30 5 

4.0 5 

6.9 5 

2.8 5 S.S 5 
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TABLE 4.17 

Comparison of experimentally observed population densities of Hel 

13 

states by Johnson with present calculations using optically thick 

conditions 


Hinnov's expt. Present cal. Present cal. par- 

values optically thigh tially optically 

thick 


n e 
n(l ) 


States 


4.0 


1.36 

8 . 2 13 


13 


4 1 S 

3.6 6 

1.8 6 

2.0 6 

5 1 S 

1. 07 6 

3.6 5 

4.7 5 

s 1 ? 

7 

1.74 / 

2.57 6 

l.l 7 

4 X P 

2. 6 6 

1.2 6 

2.1 6 

s 1 ? 

6. 6 5 

2. 4 5 

2.11 6 

3 1 D 

l.l 7 

2. 5 6 

3.4 6 

4 1 D 

2. 6 6 

1.5 5 

1.37 6 

5 1 D 

5.9 5 

1 . 5 5 

2. 0 5 

3 3 S 

1.96 7 

1.38 7 

1.57 7 

4 3 S 

4.0 6 

3.17 6 

3 . 4 6 

5 3 S 

l.l 5 

3.6 5 

4.4 5 

3 3 P 

1. 49 7 

1.6 7 

7 

2.0 

4 3 P 

2.4 6 

1.8 6 

1. 87 6 

5 3 P 

5 

■7.5 

2.8 5 

2.9 5 

3 3 D 

1.4 7 

8. 42 6 

6 

8.98° 

4 3 D 

2.6 

2.0 6 

2. 5 5 

5 3 D 

5 

6.2 

S.O 5 

5.5 5 
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Tables 4.16 and 4.17 reveal that the matching of the experi- 
mental absolute population densities with those from the 
model used in the present work is quite satisfactory. 

The present model using simultaneous solution for Hel 

and Hell levels can also be applied to the recent experiments 

in ' TPD* (Test Plasma by Direct-current Discharge) machine at 

Hagoya University, reported by Otsuka, Ikee and Ishii. For 

3 1 

Hel, absolute populations of 3 D and 4 D were measured at 

several points along the plasma column where observation 

windows are located. At these windows measurements were 

carried out also on plasma parameters n e , T £ and neutral 

density n . They did not make any attempt to explain the 
o 

observed population densities of Hel states from any theore- 
tical model. Application of the collisional— radiative model 
used in the present work to TPD machine conditions leads to 
absolute populations of all the Hel and Hell states out of 
which only the populations of 3 D and 4 D are presented in 
Table 4.18, the states for which detailed experimental data 
is available in Ref. 26. It may be noted that the present 
model provides satisfactory order of magnitude explanation of 
the experimentally observed populations. Here also, the 
partially optically thick condition with reduced rate coeffi- 
cients for the transitions 2 1 P-3 1 D are used, but at high electron 
densities, as is the case in TPD experiments, such changes in 
the rate coefficients have insignificant effect. 
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TABLE 4.18 


Comparison of experimental values of Hel 3 3 D and 4 3 D and Hell 
(1=43 population densities (Otsuka. Ikee, ishii (Ref. 263 3 with 

present calculations 

Hel 3 3 D, 4- 4 D Populations 


State T 


n 


n. 


n 


n(p) n(p) present 
expt. calculations 
Ref. 26 Opt. Opt. 

thin thick** 


3"D window 3 2.68 14 3.0 14 2.68 14 2.0 7 


3.17 eV 


1 

.4 D 


2.1 


1.5 7 2.26 7 


1 . 6 6 1 . 8 6 


14 - ^4 

p. 129 In Table 4.18, in place of n 1 values 3.0 ,1.0 , 

7.0^, read 2. 3^, 2.1^, 2.0 4 . 
p. 130, line 5. In place of heretofore, read henceforth. 


"Hell ( 4) Populations 


Window 

No. 

T 

0 

n + (4) 

observed 

n + (4) 

optically thin 

n + (4) 

thick 

Lyman 

optically 
for Hell 
radiation* 

Window 


3.17 

2.2 5 

9. 62 4 

5 

2.0“ 


Window 

4 

0.88 

2.5 6 

2. 33 7 

- 


Window 

5 

0.34 

2. 6 4 

7 

1.37 ' 

_ 



A 12 = O.OOOl, A- 13 = 0.06, A 14 = 0.1, A 15 =0.5 and A ± . = 1 


** 


k l,5 = °* 0013 ' A i,n " °*° 10 ' A i #1 8 = °* 07 ' ^,25 = °* 2 ' 


*1, 26“ °* 5 



130 


Absolute values of Hell populations were also obtained 

in the TPD machine experiments for i=4, from experimental line 

intensity measurements. The population density of this level, 

« 

in terms of a collisional-radiative model of only Hell states, 

can be written as (eg. 10’ in Ref. 26; .heretofore terminology of 

Ref. 26 is used, different from those used in the present work), 

E 1. 

n . = r (4) n^ n ++ A(4) + r, (4) n., exp C— ~ - l)] (4.42) 

4 o e 1 1 c ‘■kT . 2 

e 4 

In the terminology of Ref. 26 which considers Hell states 
separately from Hel levels but along with He ions, r Q in 
eg. (4.42) is the contribution from the He ++ continuun, r^ is. 
the contribution from the Hell ground state, n^ the population 
of Hell ground state, E^ is the ionization potential of the 
ground state of Hell, and A is given by 


A(i) = 


. 2, 3 E 

* —■ — o A> exp ( „ - ) 

(2*mkT ) J/Z i kT 

e e 


.2 4.14 x 10 

T (°K) 3//2 
e 


-16 


exp - „ - ' 
i 2 kT e 

(4.43) 


2 

In their paper, the authors miss the factor i (which has a 

O 

value of 4 =16) in the second eguaiity of eg. (4.43) and 

S 3 

calculate a value of n^=l*6xl0 cm" . Further, comparing with 

& 

the experimental value of 2.5 x 10 at window 4, they erroneously 
conclude that there is agreement between the experimentally, 
observed and theoretical values. 

Necessary correction, however, gives a theoretical 
7 

value of 2.33x10 which is off from the experimental value of 
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2.5x10 by an order of magnitude. Using appropriate r Q values 
calculated for the model used by Otsuka and others and n"*~*" 
from their Table as they have used, the values for the helium 
ion Hell (i=4) population density for the TPD plasma condi- 
tions are shown in Table 4.18. Use of Drawin's^ or McWhirter 
and Hearn's 0 Table for the Hell r^ coefficients yield essen- 
tially the same values. Clearly, the calculated n^_ values 
using Hell model chosen by Otsuka and others show large 
deviations from the experimentally observed populations. 

The conditions of their experiment however merit analysis 
in terms of applicability of optically thick conditions for 
Hell radiation. Using the same considerations as presented 
earlier, a partially optically thick model is justifiable with 
Hell Lyman radiation, optical escape factors as A^ 2 = 0.0007, 
2=0. 06, A^=0. 1, an< ^ rema i n i n 9 *-i* Applica- 

tion of r , r. values thus obtained to eg. (4.42) shows that 
for the plasma conditions of window 3, the major contribution 
comes from the ground state of Hell (i=l), hence from the r^ 
factors, whereas for the conditions at windows 4 and 5, i.e., 
at lower T g values, the dominant contribution is from the He 

continuum, i.e. from the r values. For the windows 4 and 

o 

5 cases, n concentration is required for calculation of the 
He (i=4) population, whereas for the case of window 3, n 
is not necessary. Assuming the helium ion ground state 
population as approximately equal to the electron density. 
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the population density of Hell (i=4) turns out to be 2.0xl0 5 , 
within -*v-20% of the experimentally observed value, which is at 
least as good as the value obtainable from calculations losing 
only continuum contributions as shown in Table 4.18. This 
shows that use of a partially optically thick model at T e =3.17 eV 
n 4 populations can be explained satisfactorily. 



Chapter V 
A REVIEW 


In the present work emission intensities of Hel and 
Hell line radiation have been measured in a helium plasma 
column subjected to a longitudinal magnetic field 0-700 gauss 
The measurements involve 16 lines of Hel and the- 4685 A line 
of Hell. The emission intensities were found to increase 
with increasing magnetic field and in many cases were found 
to go through a maximum and then decrease in the magnetic 
field range of the investigation. It was observed that, less 
the electron density, more pronounced is the effect of emis- 
sion intensity enhancement. 

Plasma parameters n , T were measured using Langmuir 
probes/.- n g also from emission intensities using a emission 
current ratio method, and T e also from a singlet- triplet 
Intensity ratio method. Electron densities were found to 
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increase by a factor of about 80 with increasing field, 
associated with a concurrent fall of electron temperature 
from about 11 eV to about 6 eV. 

As an attempt to explain the observed plasma line 
emission, a general collisional- radiative model is formulated 
which features completely separated sublevels up to n=5 and 
coalesced level for n=6 to n=12. 13 levels of Hell (n=l-13) 

are included in the model, mating a total of 45 levels. In 
addition to the usual radiative and electron impact processes, 
direct ionization-excitation processes to populate the excited 
states of Hell from Hel states, as well as atom-atom collisions 
are also incorporated. Neglecting diffusive losses for the 
excited states, 43 equations of Hel and Hell are simultaneously 
solved with respect to the ground state populations of Hel 
and Hell. This gives the population coefficients from which 
population densities of various levels are obtained for a 
wide range of n Q , T^ values and emission intensity enhancement 
factors are calculated therefrom. The model explains reason- 
ably well the observed enhancement factors. 

The model is then compared with collisional -radiative 
models used by other workers. It was found that the models 
which use combined singlet- triplet levels make poor predic- 
tions of enhancement factors. This is probably due to the 
marked difference in behaviour of the singlet and triplet 

levels as a function of n , T , and combination of the levels 

0 0 
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However, that the population densities which arise out of 
calculations of the collisional-radiative model (and which 
explain the emission enhancements well) is reasonable is shown 
by the extent to which it predicts the absolute population 
densities of Stellarator experiments of Johnson and TPD experi- 
ments of Otsuka, Ikee and Ishii. Experimentally measured 
population density ratios of He I and Hell states is an important 
additional point in this regard. 

As regards diagnostic methods, it is well recognized 
that no single method can be depended on for all conditions, 
in plasmas, and it is for this reason more than one method was 
used. The n (B) profile from the probe shows that it is lower 
by approximately a constant factor from the corresponding n e 
values obtained from the emission current ratio method. It is 
quite possible that a constant factor in probe density calcula- 
tion is missing - may be an estimate of the r effective .probe area, 
which if incorporated will make the probe densities close to 
the actual situations. Whereas T^ measurements from higher 
level intensities would have been desirable, the measuranents 
could not be made due to weak intensities. That weak intensi- 
ties may be expected follows from -a comparison of the population 
densities in the present work and sasme cases where such T g 
measurements have been made. In the present work, n(8)=9.93 , 

£T ET 

n(9)=5.9 , n(lO)=4.4 , whereas in Drawin, Klan and Ringler's 
work n (8) =1. 63^, n(9)=1.8^, n(10)=l,9^ and in Ikee and Takeyaraa's 
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work n(8)=8.19 8 , n(9>=s.s 8 and n(1 0)-6 21 8 re „ r k ^ 

u ; ~ b * 21 • I t ftnast be admitted 

however that a detailed seatiai n m 

x V T e prof jJL.e. Would be valuer-. 

ble. However, with present experimental arrangement very precise 
measuranents of such profiles are difficult. A 9 cm diameter 
new plasma tube with such facilities is under construction, 
where such measurements will be possible and additional side-on 
suronents of radiation, followed bv appropriate Abel trans- 
formation, ought to give spectroscopic values of local plasma 
parameters a s well. 

Whereas there are undoubted uncertainties about the 

optical escape factors, chiefly arising from the uncertainties 

of absorption lengths L, it has been clearly shown that even 

with considerable changes in the value of L, there would not be 

much difference in the matching of most of the emission enhance- 
merit factors. 

irect ionization-excitation processes included in 
the model comprise a significant addition, not known to have 
used in any other collisional -radiative model. Though 
the values of rate coefficients used are higher than the esti- 
mates of Dalgamo and McDowell, they are compatible with the 
ionization-excitation cross-section given by Drawin. It stauld 
be emphasized here that if such ionisation-excitation processes 
are not incorporated, i.e., there is no direct contribution of 
the neutral states to the ionic excited states, the observed 
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population densities of Hell cannot be explained. This is 
regardless of whether one assunes near complete ionization of ' 
He + (so that Hell states are populated mainly from the second 
^continuum and cascading) or all ions are in He + states so that 
the Hell states are populated mainly from excitation from 
Hell (i=l). if one considers that hydrogenic Hell rate coef- 
ficients are among the most accepted of the rate coef ficients, 
this point becomes more meaningful. As for the P- D rate * 
coefficient, it may be mentioned that the change is’ not too 
unreasonable compared to the degree of certainty to which most 
other rate coefficients are known. Such changes are sometimes 
made to see what makes the theoretical results correspond the 
experimental data. Furthermore, even if these coefficients 
were not changed, the basic pattern and quality of matching 
would not have been significantly worse. 

On the positive side, one might mention the following 
points about the present work. It provides a collisional- 
radiative model in which Hel and Hell levels are considered 
together, and useful applications made thereof. In higher 
temperature plasnas such model would be increasingly necessary 
where one needs to deal with more of Hell population densities. 
That it explains several experimental results of other worker's 
is an important point in support. Lastly, whereas there have 
been stray reports on emission enhancements, this is perhaps 
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the first effort, to the knowledge of the present author, 
which makes a systematic study of emission enhancements of 
helium plasmas in magnetic field, and also explains the 
observed phenomena in terms of fundamental collisional and 
radiative processes. 
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